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This article reviews recent theoretical and experimental work on transport due to the surface 
states of three-dimensional topological insulators. The theoretical focus is on longitudinal trans- 
port in the presence of an electric field, including Boltzmann transport, quantum corrections and 
weak localization, as well as longitudinal and Hall transport in the presence of both electric and 
magnetic fields and/or magnetizations. Special attention is paid to transport at finite doping, to 
the 7r-Berry phase, which leads to the absence of backscattering, Klein tunneling and half-quantized 
Hall response. Signatures of surface states in ordinary transport and magnetotransport are clearly 
identified. The review also covers transport experiments of the past years, reviewing the initial 
obscuring of surface transport by bulk transport, and the way transport due to the surface states 
has increasingly been identified experimentally. Current and likely future experimental challenges 
are given prominence and the current status of the field is assessed. 



I. INTRODUCTION 

The conventional band theory picture divides solids 
into metals, semiconductors and insulators based on the 
size of their band gaps and occupation of the conduction 
band. Within this picture, an insulator is understood 
to be a material with a band gap exceeding room tem- 
perature by orders of magnitude. In nature, other types 
of insulators exist, for example induced by a large dis- 
order concentration, as in the Anderson insulator, or by 
strong Coulomb repulsion between electrons, as in the 
Mott insulator. In this context, the discovery of topolog- 
ical insulators (TI), first posited in the groundbreaking 
research of Kane and MelejJ has been a milestone gal- 
vanizing physics researchPMS These novel materials are 
band insulators in the bulk while conducting along the 
surfaces, possessing gapless surface states with a well- 
defined spin texture protected by time-reversal symme- 
try. The landmark discovery of topological insulators has 
revolutionized the understanding of insulating behavior. 
In addition to a variety of exotic effects, Tl-based struc- 
tures may help solve long standing puzzles in physics as 
well as contribute to the development of quantum com- 
putation. Specifically, the interface between a topologi- 
cal insulator and a superconductor was predicted to sup- 
port Majorana fermionspthat is, a species of Bogolyubov 
quasiparticles represented by real fermions that are their 
own antiparticles. Topological insulators are being ex- 
plored with a view towards applications, as a potential 
platform for quantum computation^ and as rich physical 
systems in their own right. 



A. Topological order 

On a more profound level, topological insulators are a 
manifestation of topological order. Unlike other exam- 
ples of topological order, such as the fractional quan- 
tum Hall effect or chiral p-wave superconductors, the 
emergence of topologically insulating behavior is a one- 



particle phenomenon, related to the existence of strong 
spin-orbit coupling. A topologically insulating state in 
2D, the quantum spin Hall insulator, was predicted in 
HgCdTe quantum wells. The topological behavior in 
this case is related to the existence of edge states.^ The 
key concept is band inversion consider: two materials, 
of which one has an inverted band structure. In the case 
demonstrated so far this was a CdTe/HgTe quantum well. 
Each state can be labeled by a Z 2 topological invariant v 
(topological quantum number.) As the well width is in- 
creased the band gap is inverted and the Z 2 invariant of 
the ground state changes. From a lay perspective, in the 
context of the quantum spin Hall effect the Z 2 invariant 
can be thought of as a topological invariant that counts 
the number of edge states. The band inversion transi- 
tion is not characterized by symmetry breaking 121SJ ^ j g 
rather a topological phase transition: the fundamental 
symmetry of the lattice is the same on either side of the 
transition, but there is a change in the topological invari- 
ant v of the ground state. Kane and Mele^ demonstrated 
that in two dimensions the invariant v can be zero or one, 
the former representing ordinary insulators (equivalent to 
the vacuum) and the latter topological insulators. The 
two band structures corresponding to v — and v = 1 
cannot be deformed into one another. The material with 
an inverted band structure possesses a pair of crossing 
edge states which realize a quantum spin Hall insulator, 
that is, they give rise to a quantum spin Hall effect (in 
fact it possesses two pairs, by the fermion doubling theo- 
rem.) In the quantum spin Hall effect the charge current 
is carried by edge states with different spin polarizations. 
Because of spin-momentum locking due to the spin-orbit 
interaction electrons traveling along the edges in differ- 
ent directions have opposite spin orientations, and the 
chiral band structure associated with the surface states 
prohibits backscattering. The quantum spin Hall insu- 
lating state was observed experimentally shortly after its 
predictiorP. For an excellent review of this effect, which 
will not be considered here, the reader is referred to the 
work of Koenig et alW. 
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Topological insulators were also predicted to exist in 
3D EalEI where analogous sets of crossing surface states 
are present. Three-dimensional topological insulators 
will be the focus of this review. In 3D there are 
four Z2 topological invariants characterizing the band 
structure,^ which are related to the time-reversal invari- 
ant momenta and time-reversal 'polarization at eight spe- 
cial points in the Brillouin zone. The four Z 2 topological 
invariants in 3D are customarily labeled vq, v\, and 
1/3 . They can be calculated easily if the system has inver- 
sion symmetry!^ The Z 2 topological invariant vq distin- 
guishes between strong and weak topological insulators: 
weak TI have vq — whereas strong TI have vq = 1. In 
a weak TI the surface states are not protected by time 
reversal and can be localized by strong disorder. A weak 
TI can also be interpreted as a 3D stack of 2D quantum 
spin-Hall insulators. On the other hand, a strong TI is 
a new class of material and cannot be directly related 
to the quantum spin-Hall insulator. Interestingly, it was 
found that the crossover from a 2D to a 3D TI occurs 
in an oscillatory manner,^ with the material alternating 
between topologically trivial and nontrivial phases as the 
layer thickness is increased. 

In an alternative picture of TI that does not rely on Z 2 
invariants, Murakami studied the physics of gap closing 
in 2D and 3D focusing on whether the gap closes as cer- 
tain parameters are tuned. Murakami also found that 3D 
inversion-asymmetric systems possess a finite crossover 
region of gapless phase between the ordinary insulator 
and topological insulator regimes. 

The surface states of 3D TI are w ell described by a 
Rashba spin-orbit Hamiltonian.EH^Ell A detailed discus- 
sion of the model Hamiltonian for TI has been given in 
Ref. I17| and a reduced Hamiltonian including surface 
and bulk states has been presented in Ref. [T5J These 
states are chiral, possessing a well-defined spin texture, 
referred to as spin-momentum locking, and have an en- 
ergy spectrum akin to a Dirac cone. Chiral Rashba sur- 
face states have a nontrivial band topology and are char- 
acterized by a 7r Berry phase, which is associated with 
Klein tunneling, and provides protection against coher- 
ent back scattering and therefore weak localization and 
Anderson localization. 

Several materials were predicted to be topological in- 
sulators i n thr ee dimensions. The first was the alloy 
Bii—.jSba, j 19 l 2 ° l followed by the tetradymite semiconduc- 
tors Bi 2 Se3, Bi 2 Te3 and Sb 2 Te3 These materials have 
a rhombohedral structure composed of quintuple lay- 
ers oriented perpendicular to the trigonal c axis. The 
covalent bonding within each quintuple layer is much 
stronger than weak van der Waals forces bonding neigh- 
boring layers. The semiconducting gap is approximately 
0.3 eV, and the TI states exists along the (111) direc- 
tion. In particular Bi 2 Se3 and Bi 2 Te3 have long been 
known from thermoelectric transport as displaying siz- 
able Peltier and Seebeck effects, and their high quality 
has ensured their place at the forefront of experimental 
attentionP Initial predictions of the existence of chiral 



surface states were confirmed by first principles studies 
of Bi 2 Se3, Bi 2 Te3, and Sb 2 Te3.'22' Heusler alloys were re- 
cently predicted to have topological surface states,^ as 
well as chalcopyrites!^ 

Unlike graphene, the Hamiltonian of topological insu- 
lators is a function of the real spin, rather than a sub- 
lattice pseudospin degree of freedom. This implies that 
spin dynamics will be qualitatively very different from 
graphene. Moreover, the twofold valley degeneracy of 
graphene is not present in topological insulators. Despite 
the apparent similarities, the study of topological insu- 
lators is thus not a simple matter of translating results 
known from graphene. Due to the dominant spin-orbit 
interaction it is also very different from ordinary spin- 
orbit coupled semiconductors. 

Strictly speaking, the first instance of topologically in- 
sulating behavior discovered in nature is the quantum 
Hall effect (QHE). There is no difference in symmetry 
between the two sides of the quantum Hall transition, 
involving localized versus extended states. Rather, the 
Landau levels in the QHE are identified by their Chcrn 
number (TKNN invariant),^ which is the integral of the 
geometrical curvature over the Brillouin zone. From a 
mathematical perspective, the Chern number is a rela- 
tive of the genus and winding number. An analogy can 
be made between the Gaussian curvature and genus on 
the one hand and the Berry curvature and Chern num- 
ber on the other. In the QHE the Chern number can be 
thought of as the topological invariant that counts the 
number of Landau levels, being therefore equivalent to 
the filling factor. Nevertheless, quantum Hall systems 
naturally break time-reversal due to the presence of a 
magnetic field. Haldane^l devised a model of the QHE 
without an overall magnetic field and therefore without 
Landau levels, in which local magnetic fields were gen- 
erated by circulating current loops but averaged to zero 
over the whole sample. This model realizes the parity 
anomaly, but to this day remains a theoretical construct. 

In summary, in order to determine whether a certain 
material is a topological insulator it is necessary to cal- 
culate a series of topological invariants, for which one 
requires the band structure of the material. The exis- 
tence of chiral surface states is protected by time reversal 
symmetry, and is thus robust against smooth perturba- 
tions invariant under time-reversal, such as non-magnetic 
disorder and electron-electron interactions. In plain par- 
lance, in the same manner that a doughnut can be de- 
formed into a coffee cup, the band structure of a material 
cannot be deformed so as to remove the crossing of a pair 
of surface states (see the excellent figure in Ref. [H) 



B. Experimental discovery of 3D TI 

Thanks to surface-sensitive probes such as scanning 
tunneling microscopes (STM) and angle-resolved photoc- 
mission (ARPES), the existence of gapless chiral surface 
states in topological insulators has been established be- 
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yond any doubt. Experimental work originally reported 
unconventional behavior in Bi2Se3 and Bi 2 Te3,^- 7 which 
was later understood to arise from their TI surface states. 
Several years later, experimental progress on TI skyrock- 
eted. Hsieh et o/j*^ were the first to observe the Dirac 
cone-like surface states of Bii-^Sb^, with x rs 0.1, us- 
ing ARPES. Following that, the same group used spin- 
resolved ARPES to measure the spin polarization of the 
surface states, demonstrating the correlation between 
spin and momentum.^ Shortly afterwards experiments 
identified the Dirac cones i n the tetradymite semicon- 
ductors Bi 2 Se 3 p2HHI Bi 2 Te 3 | ™> and Sb 2 Te 3 P Reflec- 
tion high-energy electron diffraction (RHEED) was used 
to monitor MBE sample growth dynamics of Bi 2 Te3,^ 
and the electronic structure has been shown to be control- 
lable in thin filmsP^TI growth has been covered in recent 
reviews On a related note, recently Bi 2 Se3 nanowires 
and nanoribbons have also been grownj^ and a topolog- 
ical phase transition was reported in BiTl(Si_dSe £ ;)2}2^' 
which may lead to the observation of charge e/2. 

Zhu et alW^ used ARPES to study the surface states 
of Bi 2 Se3. These researchers employed in-situ potassium 
deposition to tackle the instability of electronic proper- 
ties, which otherwise evolve continuously in time. In their 
samples, as potassium is deposited, new spin-polarized 
quantum- well states emerge (with no k z dispersion), orig- 
inating in the parabolic bulk bands and distinct from 
the Rashba-Dirac cone. These new states also follow a 
Rashba-like dispersion, and their spin splitting is tunable 
as well as reversible. The potassium-induced potential 
gradient dV/dz enhances the spin splitting. These states 
may affect transport data even from pristine surfaces. 

Considerable efforts have been devoted to scanning 
tunneling microscopy (STM) and spectroscopy (STS), 
which enable the study of quasiparticle scattering. Scat- 
tering off surface defects, in which the initial state inter- 
feres with the final scattered state, results in a standing- 
wave interference pattern with a spatial modulation de- 
termined by the momentum transfer during scattering. 
These manifest themselves as oscillations of the local den- 
sity of states in real space, which have been seen in several 
materials with topologically protected surface states. 

Roushan et al™ imaged the surface states of 
Bi .g 2 Sb .o8 using STM and spin-resolved ARPES 
demonstrating that scattering between states of opposite 
momenta is strongly suppressed in the presence of time- 
reversal invariant disorder. Gomes et alW^ investigated 
the (111) surface of Sb, which displays a topological metal 
phase, and found a similar suppression of backscattering. 
Zhang et alW^ and Alpichshev et aiP^ observed the sup- 
pression of backscattering in Bi 2 Te3, and used STM to 
image surface bound states.^ The absence of backscat- 
tering has been studied theoretically in the conte xt of 
quasiparticle interference seen in STM experiments! 45 * 46 ^ 
Park et alW^ performed ARPES accompanied by theo- 
retical modelling on Bi 2 Se3. The observed quasiparticle 
lifetime of the surface states is attributed to quasiparticle 
decay predominantly into bulk electronic states through 



electron-electron interaction and defect scattering. Stud- 
ies on aged surfaces show the surface states to be unaf- 
fected by adsorbed atoms or molecules on the surface, 
indicating protection against weak perturbations. 

Experimental studies have provided evidence of the ex- 
istence of chiral surface states and of their protection 
by time-reversal symmetry. Experimental research has 
begun to investigate the fascinating question concerning 
the transition between the three- and two-dimensional 
phases in this material.^ Zhang et alW^ used ARPES to 
study a thin 3D Bi 2 Se3 slab grown by MBE, in which 
tunnelling between opposite surfaces opens a small, 
thickness-dependent gap. The energy gap opening is 
clearly seen when the thickness is below six quintuple 
layers. The gapped surface states also exhibit sizeable 
Rashba spin-orbit splitting due to the potential differ- 
ence between the two surfaces induced by the substrate. 
In other words, one of the most important findings of 
this experiment is that thin films need to be thicker than 
approximately six quintuple layers, or 6nm. Bianchi et 
aZ.HS showed that the band bending near the Bi 2 Se3 sur- 
face leads to the formation of a two-dimensional electron 
gas (2DEG), which coexists with the topological surface 
state in Bi 2 Se3. In the setup used by this experimental 
group, a topological and a non-topological metallic state 
are confined to the same region of space. 

Several efforts have focused on the role of magnetic im- 
purities on the surface states of 3D TI. The work of Hor et 
al. demonstrated that doping_Bi 2 Te3 with Mn results in 
the onset of ferromagnetism. 50 Chen et al. 51 introduced 
magnetic dopants into Bi 2 Se3 to break the time reversal 
symmetry and further position the Fermi energy inside 
the gaps by simultaneous magnetic and charge doping. 
The resulting insulating massive Dirac fermion state was 
observed by ARPES. More recently, Wray et oi.^ studied 
the topological surface states of Bi 2 Se 3 under Coulomb 
and magnetic perturbations. The authors deposited Fe 
on the surface of Bi 2 Se3, using the fact that Fe has a 
large Coulomb charge and sizable magnetic moment to 
modify the spin structure of the Bi 2 Se3 surface states. 
Interestingly, it was found that this perturbation leads 
to the creation of odd multiples of Dirac fermions. In 
other words, new pairs of Dirac cones appear as more Fe 
is deposited, also described by Rashba interactions, and 
the interaction grows in strength as Fe is added. Numeri- 
cal simulations indicate that the new Dirac cones are due 
to electrons localized deeper inside the material. At the 
same time a gap is opened at the original Dirac point, as 
one expects. 



C. The importance of transport 

The key to the eventual success of topological insu- 
lators in becoming technological materials is inherently 
linked to their transport properties. Potential applica- 
tions of topological surface states necessarily rely on the 
realization of an edge metal (semimetal), allowing con- 
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tinuous tuning of the Fermi energy through the Dirac 
point, the presence of a minimum conductivity (maxi- 
mum resistivity) at zero carrier density, and ambipolar 
transport. At this point in time experimental and theo- 
retical studies of equilibrium TI abound. Transport has 
been lagging somewhat but is now picking up in both 
theory and experiment at a brisk rate. 

Despite the success of photoemission and scanning tun- 
neling spectroscopy in identifying chiral surface states, 
signatures of the surface Dirac cone have not yet been 
observed in transport. Simply put, the band structure 
of topological insulators can be visualized as a band in- 
sulator with a Dirac cone within the bulk gap, and to 
access this cone one needs to ensure the chemical poten- 
tial lies below the bottom of the bulk conduction band. 
Given that the static dielectric constants of materials un- 
der investigation are extremely large, approximately 100 
in Bi2Se3 and 200 in Bi2Te3, gating in order to bring 
the chemical potential down is challenging. Therefore, in 
all materials studied to date, residual conduction from 
the bulk exists due to unintentional doping. The above 
presentation makes it plain that, despite the presence 
of the Dirac cone, currently no experimental group has 
produced a true topological insulator. Consequently, all 
existing topological insulator systems are in practice ei- 
ther heavily bulk-doped materials or thin films of a few 
monolayers (thus, by definition, not 3D.) Following the 
pattern set by the imaging of the surface Dirac cone, the 
Landau levels ob served in STM studies have not been 
seen in transport! 53 * 54 ! These last works illustrate both 
the enormous potential of the field and the challenges to 
be overcome experimentally, on which more below. 

In this work I will provide a theoretical framework for 
understanding TI transport, review existing TI transport 
theories, and survey the impressive experimental break- 
throughs in transport registered in the past years. The 
review will focus on phenomena in electric fields as well 
as electric and magnetic fields, but not magnetic fields 
alone. As it stands, fascinating aspects such as the search 
for Majorana fermions also lie outside the scope of this re- 
view, since no transport work exists on these topics. The 
transport work covered consists of two main parts, which 
in turn can be subdivided into five smaller parts: or- 
dinary (non-magnetic) transport, comprising Boltzmann 
transport and weak localization, and magnetotransport, 
comprising the quantum Hall, ordinary Hall and anoma- 
lous Hall effects. Theoretical work in these fields will 
be covered first, followed by experiment, and I will close 
with a section on experimental challenges, followed by 
general conclusions and an outline of future directions. 



II. EFFECTIVE HAMILTONIAN 

The effective Hamiltonian describing the surface states 
of non-magnetic topological insulators can be written in 



the form of a Rashba spin-orbit interaction^^!] 

H ok = Acr ■ k x z = -Ak cr ■ 9. (1) 

It is understood that k = (k x ,k y ), and 9 is the tan- 
gential unit vector corresponding to k. The constant A 
= 4.1eVl for Bi 2 Se 3 P The eigenenergies are denoted 
by e± = ±Ak, so that the spectrum possesses particle- 
hole symmetry. The Hamiltonian is similar to that of 
graphene, with the exception that the vector of Pauli 
spin matrices cr represents the true electron spin. A 
scalar term Dk 2 also exists in principle. The scalar and 
spin-dependent terms in the Hamiltonian are of the same 
magnitude when Dkp = A, corresponding to kp ss 10° 
m _1 in Bi2Se3, that is, a density of 10 14 cm -2 , which 
is higher than the realistic densities in transport exper- 
iments. Consequently, in this review scalar terms will 
not be considered. I will also not take into account small 
anisotropy terms in the scalar part of the Hamiltonian, 
such as those discussed in Ref. ETJ Terms cubic in k in 
the spin-orbit interaction are in principle also present, in 
turn making the energy dispersion anisotropic, but they 
are much smaller than the fc-linear terms in Eq. ([!]) and 
are not expected to contribute significantly to quantities 
discussed in this review. All theories described here re- 
quire Eft/H 3> 1 to be applicable, with r the momentum 
relaxation time. 

A topological insulator can doped with magnetic im- 
purities, which give a net magnetic moment as well as 
spin-dependent scattering. Magnetic impurities are de- 
scribed by an additional interaction Hamiltonian 

^mag(r) =a-Y l V(r-R I )a I , (2) 
/ 

where the sum runs over the positions Rj of the mag- 
netic ions with sj. To capture the physics discussed 
in this review it is sufficient to approximate the poten- 
tial V(r — Rj) = J 5(r — Ri), with J the exchange 
constant between the localized moments and itinerant 
carriers. The magnetic ions are assumed spin-polarized 
so that sj = sz. Fourier transforming to the crystal 
momentum representation {|fc)}, the fc-diagonal term, 
i?,^g = n mag J s<T z = M u Z: gives the magnetization, 
with n mag the density of magnetic ions. The fc-off- 
diagonal term causes spin-dependent scattering 

= y« z ^ {k ~ k ' yRl > (3) 

with V the volume. The effective band Hamiltonian in- 
corporating a mean-field magnetization takes the form 

H mk = -Akcr-9 + crM=^cr-fl k . (4) 

Retaining the same notation as in the non-magnetic case, 
the eigenenergies are e± = ±V A 2 k 2 + M 2 . The bulk of 
the work on magnetic TI has concentrated on out-of-pane 
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magnetizations and magnetic fields. An in-plane mag- 
netization or magnetic field does not alter the Rashba- 
Dirac cone, though coherence between layers can induce 
a quantum phase transition.^ 

Interaction with a static, uniform electric field E is 
quite generally contained in H E ,kk> = H Ekk , + H S J kk ,, 
the scalar part arising from the ordinary position op- 
erator H s E c kk , — (eE ■ r)fcfc'l, with 1 the identity ma- 
trix in spin space, and the side-jump part H^ kk , = 
eXcr ■ (fe x E) 5 kk i arising from the spin-orbit modifica- 
tion to the position operator, r — > f + Act x fej^ with A 
a material-specific constant. 

Elastic scattering off charged impurities, static defects 
and magnetic ions, but not phonons or other electrons, 
is contained in the disorder scattering potential Ukk f = 
Ukk'Yl j e > ( k ~ k ')' Rj with Rj the random locations of 
the impurities. The potential Ukk' incorporates the mag- 
netic scattering term 

H^g ofEq. H). The potential of 
a single impurity Ukk' = (I — iAer-fcx k')lAkk' — Jsa z , 
with Ukk' the matrix element of a Coulomb potential be- 
tween plane waves, and Xkp <C 1. Note that Ukk' has 
the dimensions of energy times volume. 

Charged impurities are described by a screened 
Coulomb potential. The screening function can be eval- 
uated in the random phase approximation (RPA.) In 
this approximation the polarization function is obtained 
by summing the lowest bubble diagram, and takes the 
forrrP (A, A' = ±) 



I%,w) = - 



./'« 



OfcA 



Ofc'A' 



A ^ u 

k\\' 



£k\ — £fc'A' 



IT) 



1 + AA' cos 7 



(5) 

where fok\ = faiskx) is the equilibrium Fermi distribu- 
tion function. The static dielectric function, which is of 
relevance to the physics discussed in this review, can be 
written as e(q) = 1 + v(q) IT(q), where v(q) — e 2 /(2eoe r g), 
where e r is the relative permittivity. To determine 11(g), 
we assume T = and use the Dirac cone approxima- 
tion as in Ref. 57, This approximation is justified in 
the regime T/Tp <C 1, with Tp the Fermi tempera- 
ture. At T — for charged impurity scattering the long- 
wavelength limit of the dielectric function ia^ 



<<?) = 1 



' kp 

4:TT€oe r A \ q 



(6) 



yielding the Thomas- Fermi wave vector as krF = 
e 2 kp/(4:TTeoe r A). As a result, in topological insulators 
the matrix element Ukk' of a screened Coulomb poten- 
tial between plane waves is given by 



Ukk' 



1 



Ze 2 



2e e r \k — k'\ + krF ' 



(7) 



where Z is the ionic charge (which I will assume for sim- 
plicity to be Z = 1) and krF is the Thomas-Fermi wave 
vector. The polarization function was also calculated in 
Ref. 1581 The Wigner-Seitz radius r s , which parametrizes 
the relative strength of the kinetic energy and electron- 
electron interactions, is a constant for the Rashba-Dirac 
Hamiltonian, and is given by r s = e 2 /(2ireoe r A). 

The full Hamiltonian H k = H ok + H E kk' + Ukk' + 
■H^iag • The current operator j has contributions from 
the band Hamiltonian, 
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(8) 



as well as from the applied electric field j E = c x 
E, and from the disorder potential ju = 2 ^ er x (k — 
k')Ukk>- These latter two cancel on physical grounds, as 
they represent the force acting on the system. 

Transport experiments on topological insulators seek 
to distinguish the conduction due to the surface states 
from that due to the bulk states. As mentioned in the in- 
troduction, it has proven challenging to lower the chemi- 
cal potential beneath the bottom of the conduction band, 
so that the materials studied at present are not strictly 
speaking insulators. It is necessary for ep to be below the 
bulk conduction band, so that one can be certain that 
there is only surface conduction, and at the same time 
I recall the requirement that spr/h 3> 1 for the kinetic 
equation formalism to be applicable. These assumptions 
will be made throughout the theoretical presentation out- 
lined in what follows. 



III. TRANSPORT THEORY FOR 
TOPOLOGICAL INSULATORS 



As stated above, work on transport in topological in- 
sulators has begun to expand energetically. The focus 
of experimental research in particular has shifted to the 
separation of bulk conduction from surface conduction, 
and, ideally, strategies for the elimination of the former. 
Hence it is necessary to know what contributions to the 
conductivity are expected from the chiral surface states. 

Transport problems are customarily approached within 
the framework of the linear response formalism. The re- 
sulting Kubo formula is frequently expressed in terms of 
the Green's functions of the system. This can be done 
in the Keldysh formulation, which has been covered in 
countless textbooks, one of the clearest and most com- 
prehensive being Ref. 59, Most generally, the finite- 
frequency conductivity tensor is given by 
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In the above formula, n stands for the carrier number 
density, Tr is the operator trace, () denotes the average 
over impurity configurations, po is the equilibrium den- 
sity matrix, v is the velocity operator, and G R ' A are the 
usual retarded and advanced Green's functions. In topo- 
logical insulators the indices a, ft = x, y. One important 
feature of the Kubo formula is that the conductivity is 
effectively given by the impurity average of the prod- 
uct of two Green's functions, which can be iterated in 
the strength of the disorder potential, or alternatively 
in powers of the quantity H/(eft), which as we recall 
is required to be small for the transport theory outlined 
here to be valid. This iteration is encompassed within 
the Bethe-Salpeter equation, and can be expressed di- 
agrammatically. Depending on the approximation re- 
quired, several sets of diagrams must be summed. In a 
non-magnetic system, the contribution to leading-order 
in K/(eft) recovers the Boltzmann conductivity, and cor- 
responds to the sum of ladder diagrams. Since the sum 
of ladder diagrams can be interpreted classically as rep- 
resenting diffusion of carriers, this terms also bears the 
name of diffuson. In the next order, which yields a term 
independent of r, the weak localization contribution to 
the conductivity is obtained as the sum of maximally- 
crossed diagrams, or the Cooperon (plus some minor cor- 
rections). Weak localization is usually unimportant for 
Eft/H ^> 1 unless a magnetic field is applied to measure 
the magnetoresistance. However when the Hamiltonian 
is spin-dependent/chiral it is possible to have terms in- 
dependent of t even in the weak momentum scattering 
limit. Finally, it must be borne in mind that, in topolog- 
ical insulators, identifying and explaining all the relevant 
contributions to transport necessitates a matrix formula- 
tion, thus in Eq. (|9| the equilibrium density matrix, the 
velocity operators and Green's functions are matrices. 

The Kubo formula is extremely well known in all its 
forms, and there appears to be little motivation for this 
review to focus on conventional linear response. An al- 
ternative matrix formulation, which contains the same 
physics and is potentially more transparent, relies on the 
quantum Liouville equation to derive a kinetic equation 
for the density matrix. This theory was first discussed for 
graphene monolayers^ and bilayersj^ and was recently 
extended to topological insulators including the full scat- 
tering term to linear order in the impurity density.^! Pe- 
culiarities of topological insulators, such as the absence 
of backscattering, which reflects the 7r Berry phase and 
leads to Klein tunneling, are built into this theory in a 
transparent fashion. 

The derivation of the kinetic equation for a system 
driven by an electric field in the presence of random, un- 



corrclatcd impurities begins with the quantum Liouville 
equation for the density operator pJ^H 



dp 



(10) 



where the full Hamiltonian H = Hji+ He + U, the band 
Hamiltonian H Q is defined in Eq. M, He — eE ■ r repre- 
sents the interaction with external fields, r is the position 
operator, and U is the impurity potential. We consider 
a set of time-independent states {|fes)}, where k indexes 
the wave vector and s the spin. The matrix elements of 
p are pkk' = Pkk' = (k s \pW s '), with the understand- 
ing that pkk' is a matrix in spin space. The terms H 
and He are diagonal in wave vector but off-diagonal in 
spin, while for elastic scattering in the first Born approx- 
imation Uj* k , — Ukk'^ss'- The impurities are assumed 
uncorrelated and the average of (ks\U\k's')(k's'\U\ks) 
over impurity configurations is (rii\Ukk'\ "ss') /V , where 
rii is the impurity density, V the crystal volume and Ukk' 
the matrix element of the potential of a single impurity. 

The density matrix pkk 1 is written as pkk' = fk &kk' + 
9kk' i where fk is diagonal in wave vector (i.e. fkk' oc 
5kk') while gkk' is off-diagonal in wave vector (i.e. k ^ k' 
always in gkk'-) The quantity of interest in determining 
the charge current is fk, since the current operator is 
diagonal in wave vector. We therefore derive an effec- 
tive equation for this quantity by first breaking down the 
quantum Liouville equation into 



dfk 
dt 

dgkk' 
dt 



- h [H 0k ; 



fk] = 



-~[ffjf,A]-~[Muftla) 



T W,g]kk' 



III 



kk' 



(lib) 



The term -- [U,g] k k on the RHS of Eq. (Hal will be- 



come the scattering term. Solving Eq. ( lib ) to first order 
in the scattering potential U 



gkk' 



dt 1 e - im '' h 



U,f(t-t') 



jHt'/h 



kk' 



We are interested in variations which are slow on the scale 
of the momentum relaxation time, consequently we do 
not take into account memory effects and f(t—t') ~ f(t). 
Time integrals such as the one appearing in Eq. ( [l2| are 
performed by inserting a regularizing factor e _,)t and 
letting rj — > in the end. It is necessary to carry out an 
average over impurity configurations, keeping the terms 
to linear order in the i mpuri ty density rij. 

The solution of Eq. (lib) to first order in U repre- 



sents the first Born approximation, which is sufficient 
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for ordinary Boltzmann transport, with or without mag- 
netic interactions. In order to recover phenomena such 
as skew scattering, elaborated upon below, one needs to 
work to second order in U. The total scattering term 
J(fk) = J Born (fk) + J 3rd (fk), with 

J Born (h) = iU, a-^[U, f] e^]) , (13) 



while () represents averaging over impurity configura- 
tions. In the absence of scalar terms in the Hamilto- 
nian the scalar and spin-dependent parts of the density 
matrix, n k and S k , are decoupled in the first Born ap- 
proximation. In the second Born approximation, which 
is needed in magnetic systems, one obtains the additional 
scattering term 



J 6rd (fk) 



h 3 



dt' 



dt"[U, e-*-sr[U, 



[UJ] 



(14) 



fefe 



The former can be further broken down into J 



Born 



(fk) 



be written as 



Mfh) + Jf s orn (h) + J Born (h)- The first term, J (/ fc ), 
in which A = E = in the time evolution opera- 
tor, represents elastic, spin-independent, pure momen- 
tum scattering. In Jf s orn (fk) we allow A to be finite 
but E = 0. In J Born (fk) both A and E are finite, thus 

Jfj° rn {fk) acts on the equilibrium density matrix f ok - 
Beyond the Born approximation we retain the leading 
term J 3rd {f k ) = J s 3 J d (/fc), with A finite but E = 0, and 
which is customarily responsible for skew scattering^. 
The contribution due to magnetic impurities is also oc er 2 
and is contained in Jf s orn {f k )- 

First consider the case in which the magnetization 
M = 0. In this case all terms oc A above are zero, and 
it is safe to take A = 0. To first order in E, the diagonal 
part fk satisfies^ 



dfk 
dt 



[H kJk] + J(fk) = V k . 



(15) 



The driving term V k = -| [H% , f ok ], and f ok is the 
equilibrium density matrix, which is diagonal in k. The 
fc-diagonal part of the density matrix f k is a 2 x 2 Hermi- 
tian matrix, which is decomposed into a scalar part and 
a spin-dependent part. One writes f k = n k 1L + S k , where 
S k is a 2 x 2 Hcrmitian matrix, which represents the 
spin-dependent part of the density matrix and is written 
purely in terms of the Pauli a matrices. Every matrix 
in this problem can be written in terms of a scalar part, 
labeled by the subscript n, and <x. Rather than choos- 
ing Cartesian coordinates to express the latter, it is more 
natural to identify two orthogonal directions in recipro- 
cal space, denoted by £l k = —9 and k. One projects a as 
<7fc II = <t Hfc and a k> ± = cr k. Note that <J k} \\ commutes 
with H ok , while <r kt ± does not. One projects S k onto the 
two directions in fc-space, writing S k = S^ai + S ki ±, and 
defines 5 fe ,|| = (1/2) s k ,\\ <r k .\\ and S kt ± = (1/2) s k ,j_ a k . X - 
There is no coupling of the scalar and spin distributions 
because of the particle-hole symmetry inherent in the 
Rashba-Dirac Hamiltonian. 

Introducing projection operators Pm and P± onto the 
scalar part, a k \\ and a kt ± respectively, equation (15) can 



dSu 



dt 



P\\J(fk 



dS, 



fe,_L 



dt 



[H k ,S kX } + P±J(fk) 



- Z? fc ,||, (16a) 
= V h ,±. (16b) 



The projector Py acts on a matrix M. as tr (M.a k n), with 
tr the matrix trace, while P± singles out the part orthog- 
onal to -ffofc. 

Let 7 represent the angle between k and k! . The fol- 
lowing projections of the scattering term acting on the 
spin-dependent part of the density matrix are needed 



P\\J(Sk\\) = 



P±J(Sk\\ 



P\\J{S k ±) 



krii (Tfe|| 
8hnA 

8hirA 

krii cr fc jj 
8hnA 



d9'\U kk >\ 2 (s kl \-s k ,n)(l + co8 7 ) 



dO' \ U kk >\ 2 (sfeii - s fe /|i) sin 7 



dO' \U kk >\ 2 (s fc _L + Sfe'x) sin 7. 

(17) 

Notice the factors of (1 + cos 7) and sin 7 which pro- 
hibit backscattering and give rise to Klein tunneling. In 
fact, the scattering integrals originally contain energy 8- 
functions S(e' ± — e±) = (1/A) S(k' — k), which have al- 
ready been integrated over above. The 5-functions of 
S- are needed in the scattering term in order to ensure 
agreement with Boltzmann transport. This fact is al- 
ready seen for a Dirac cone dispersion in graphene where 
the expression for the conductivity found in Ref. 1601 using 
the density-matrix formalism agrees with the Boltzmann 
transport formula of Ref. 63 (the definition of r differs by 
a factor of two in these two references.) The necessity of 
keeping the e_ terms is a result of Zitterbewegung: the 
two branches e± are mixed, thus scattering of an electron 
requires conservation of e+ as well as e_ . 

For M = the equilibrium density matrix pofc = 
(1/2) (/ ofe + + /on-) 1 - (1/2) (f ok + - fok-) a- ■ 0. One 
decomposes the driving term into a scalar part T> n (not 
given) and a spin-dependent part T> s . The latter is fur- 
ther decomposed into a part parallel to the Hamiltonian 
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T>» and a part perpendicular to it T>± 



1 eE 

2 ~li 

1 eE- d 



9/, 



Ofe- 



lo-fcll 



2 1 



a k\\ 



2 ftfc 



1 



ofe -i 



fok-) (J k± — 77 dk± <Jk±- 



(18) 

Equation ( 16 ) must be solved perturbatively in the small 



parameter WJ(eft) (which is proportional to the impu- 
rity density rii) as described in Ref. In transport 
this expansion starts at order [h/ (eft)]^ 1 (i.e. n~ l ), re- 
flecting the competition between the driving electric field 
and impurity scattering resulting in a shift of the Fermi 
surface. The leading terms in [h/ '(eft)]" 1 appear in n^ 
and 5fc iiP^l To leading order in /i/(£ft) one finds 



S 



Ek, 



dk 



1 

T 



reE-k <9/ ( 

Ah 
krii 



■ a fc|| 



27 



(19) 



|i7" fcfc /| 2 sin 2 7. 



The sin 2 7 represents the product (1 + cos7)(l — COS7). 
The first term in this product is characteristic of TI and 
ensures backscattering is suppressed, while the second 
term is characteristic of transport, eliminating the effect 
of small-angle scattering. The conductivity linear in r 
arises from the parallel part of the density matrix and is 
given by 



Boltz 



e 2 Ak F T 
~h Ah ' 



(20) 



We emphas ize ag ain that this result tends to the result 
for graphen d 60 * 63 * when D — > (note that the definition of 
t in Rcfs. 60 63 differs by a factor of 2.) The contribution 
to the conductivity independent of r is is zero for finite 
/x, which is the case considered here. Thus, the leading 
order term in the conductivity is cx r, and in the limit 
£pT jh ^> 1 there is no term of order (r)°. 

The charge current is proportional to the spin operator, 
as can be seen from Eq. ^ . Therefore a nonzero steady- 
state surface charge current automatically translates into 
a nonzero steady-state surface spin density. The spin 
density response function can easily be found by simply 
multiplying the charge current by h 2 / (—2eA) , yielding a 
spin density of 



eE x AkpT 

An 4 



(21) 



For a sample in which the impurities are located on the 
surface, with n/n, = 0.5 and E x = 25000V/m, the spin 
density is ~ 5 x 10 14 spins /m 2 (where spin = H), which 
corresponds to approximately 10 -4 spins per unit cell 
area. This number, although small, can be detected ex- 
perimentally using a surface spin probe such as Kerr ro- 
tation. It is also a conservative estimate: for very clean 
samples, having either a smaller impurity density or im- 
purities located further away from the surface, this num- 
ber can reach much higher values. Current-induced spin 



polarization is a definitive signature of two-dimensional 
surface transport in topological insulators since there is 
no spin polarization from the bulk: bulk spin densities 
vanish if the bulk has inversion symmetry, which is the 
case in the materials Bi2Se3 and Bi2Te3pH 

I will discuss next implications of these results for 
transport experiments!^! Below I will use the terminol- 
ogy high density as meaning that the density is high 
enough that a scalar term in the Hamiltonian, of the form 
Dk 2 , has a noticeable effect (meaning n » 10 13 cm -2 and 
higher) , with low density reserved for situations in which 
this term is negligible. The conductivity is a function of 
two main parameters accessible experimentally: the car- 
rier number density and the impurity density /scattering 
time. The dependence of the conductivity on the carrier 
number density arises through its direct dependence on 
kF and through its dependence on n through r. In terms 
of the number density the Fermi wave vector is given by 
kp = Ann. The number-density dependence of r depends 
on the dominant form of scattering and whether the num- 
ber density is high or low as defined above. For charged 
impurity scattering the conductivity contains terms oc n 
and n 3 / 2 (the latter being due to the scalar term not given 
in detail here.) At the same time, in a two-dimensional 
system surface roughness gives rise to short-range scat- 
tering as discussed in Ref. |M1 For short range scattering 
the two terms in the conductivity are a constant and n 1 / 2 
(the latter again due to the scalar term.) These results 
are summarized in Table |TJ where I have listed only the 
number density dependence explicitly, replacing the con- 
stants of proportionality by generic constants. 

The frequently invoked topological protection of TI 
surface transport represents protection only against lo- 
calization by back scattering, not protection against im- 
purity or defect scattering, i.e. against resistive scatter- 
ing in general. The presence of impurities and defects 
will certainly lead to scattering of the TI surface carriers 
and the surface 2D conductivity will be strongly affected 
by such scattering. If such scattering is strong (which 
will be a non-universal feature of the sample quality), 
then the actual surface resistivity will be very high and 
the associated mobility very low. There is no guaran- 
tee or protection against low carrier mobility in the TI 
surface states whatsoever, and unless one has very clean 
surfaces, there is little hope of studying surface transport 
in 2D topological states, notwithstanding their observa- 
tion in beautiful band structure measurements through 
ARPES or STM experiments. Strong scattering will lead 
to low mobility, and the mobility is not a topologically 
protected quantity. To see surface transport one there- 
fore requires very clean samples. 

At carrier densities such that Eft/H 3> 1 the conduc- 
tivity is given by <r x ° ltz from Eq. (20 1. As the number- 



density is tuned continuously through the Dirac point, 
ri — > and the condition eft/Ii 3> 1 is inevitably vio- 
lated, and a strong renormalization of the conductivity 
is expected due to charged impurities JSSHSZ! Quite gener- 
ally, charged impurities give rise to an inhomogeneous 
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TABLE I: Carrier density dependence of the longitudinal 
conductivity a xx . The numbers a\ — 0,4 represent constants. 
The table is adapted from Ref. 1621 
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FIG. 1: Calculated conductivity limited by screened charged 
impurities as a function of the surface carrier density for dif- 
ferent impurity locations d = 0, 5, 10A. In this figure we 
use the following parameters: an impurity density rn = 10 13 
cm -2 , A = 4.1 eVA, which corresponds to the Fermi veloc- 
ity vp = 6.2 x 10 7 cm/s, and the static dielectric constant 
e r — 100. The figure is adapted from Ref. 1621 



Coulomb potential, which is screened by both electrons 
and holes. The net effects of this potential are an in- 
homogeneity in the carrier density itself and a shift in 
the Dirac point as a function of position. At high den- 
sities the spatial fluctuations in the carrier density are 
of secondary importance, and do not modify the linear 
dependence of the conductivity on n, yet as the chemical 
potential approaches the Dirac point, where the average 
carrier density (n) = 0, these fluctuations play the dom- 
inant role in conduction.^ At low densities, the carriers 
cluster into puddles of electrons and holes, and a residual 
density of carriers is always present, making it impossible 
to reach the Dirac point experimentally.^ Consequently 
the renormalization due to the presence of electro n and 
hole puddles displays a strong sample dependence! 63 * 66 * 
In addition, exchange and correlation effects make a sig- 
nificant contribution to this minimum conductivity.^ An 
accurate determination of this enhancement for topolog- 
ical insulators requires detailed knowledge of the impu- 
rity density distribution, which can only be determined 
experimentally or modeled numerically by means of e.g. 



an effective medium theory.^ However, a self consistent 
transport theory provides a physically transparent way 
to identify the approximate minimum conductivity. The 
magnitude of voltage fluctuations can be calculated in 
the random phase approximation, and the result used to 
determine the residual density and the critical number 
density at which the transition occurs to the regime of 
electron and hole puddles, where the carrier density will 
be highly inhomogeneous. Both the residual density and 
the critical density are proportional to rii by a factor of 
order unityj^ and as a first approximation one may take 
both of them as « n^. In view of these observations, 
the minimum conductivity plateau will be seen approx- 
imately at <Jx° ltz , in which with the carrier density n is 
replaced by n^. Using the expressions found above for 
<Jx° ltz and t, this yields for topological insulators 



(22) 



where the dimensionless It c is determined by r s .E21With 
r s much smaller than in graphene due to the large dielec- 
tric constant, the minimum conductivity may be substan- 
tially larger!^ (Note that for a Rashba-Dirac cone the 
definition of r s contains some arbitrariness!^!) 

Transport measurements require the addition of metal- 
lic contacts on the surface of the topological insulator. 
Since the properties of the surface states depend cru- 
cially on the boundary conditions, the natural question 
is how these properties are affected by the metallic con- 
tacts. To address this question numerical calculations 
for the minimal tight-binding model studied in Ref. [55] 
were performed, focusing on the case of the case of large 
area contacts.^ When the TI surface is in contact with a 
metal, the surface states penetrate inside the metal and 
hybridize with the metallic states. A typical hybridized 
state is shown in Fig. [2^,. The amplitude of this state 
near the boundary of the topological insulator is reduced 
by a factor 1/Lg with respect to the amplitude of a 

pure surface state with the same energy, where L z is 
the width of the metal in the direction perpendicular to 
the interface. However, the local density of states near 
the boundary is not reduced, as the number of hybridized 

states also scales with . The spectrum of a topologi- 
cal insulator in contact with a metal is shown in Fig. [2]d. 
For each state the total amplitude within a thin layer of 
topological insulator in the vicinity of the boundary was 
calculated, shown in Fig. ^jp. Instead of the sharply de- 
fined Dirac cone that characterizes the free surface, one 
has a diffuse distribution of states with boundary contri- 
butions. The properties of the interface amplitude distri- 
bution shown in Fig. |2j;, i.e., its width and the dispersion 
of its maximum, depend on the strength of the coupling 
between the metal and the topological insulator. If the 
distribution is sharp enough and the dispersion does not 
deviate significantly from the Dirac cone, the transport 
analysis presented above should be applied using bare pa- 
rameters (i.e. parameters characterizing the spectrum of 
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FIG. 2: (a) Amplitude of a metallic state hybridized with a 
topological insulator surface state as function of distance from 
the interface (in units of interlayer spacing), (b) Spectrum 
of a topological insulator in contact with a metal (black). 
The states of a topological insulator with a free surface are 
also shown (purple - bulk states, red - surface states), (c) 
Density plot of the total amplitude within a thin region of 
topological insulator in the vicinity of the interface (red/dark 
grey area in panel a). Note that, instead of a well-defined 
surface mode, one has a diffuse distribution of states with 
boundary contributions. The dispersion of the maxima of 
this distribution is represented by green points in panel (b). 
The figure is adapted from Ref. 62. 



a free surface.) Nonetheless, significant deviations from 
the free surface dispersion are possible. For example, the 
location of the Dirac point in the model calculation con- 
sidered here is fixed by symmetry. However, in real sys- 
tems the energy of the Dirac point can be easily modified 
by changing the boundary conditions. Consequently, the 
effective parameters (e.g. A — »• A c g) entering transport 
coefficients may differ significantly from the correspond- 
ing parameters extracted from ARPES measurements. 
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FIG. 3: Contributions to the Hall conductivity of a TI thin 
film due to the orbital magnetic moment plotted as a function 
of Fermi level. The two plots are for film thicknesses of (a) 
20A and (b) 32A respectively. This effect can be understood 
as a spin-Hall effect in much the same manner as the 2D quan- 
tum spin-Hall insulator, and is accompanied by a topological 
phase transition as the film thickness is changed. This figure 
has been adapted from Ref. 1691 



In topological insulator thin films tunneling may be 
possible between the top and bottom surfaces. Evidently 
inclusion of this physics requires one to start with a 4 x 4 
Hamiltonian, rather than the 2x2 model discussed so 
far. Interlayer hybridization induces a gap in the surface 
states of thin films, as was discussed in Ref. [Ml Once 
the gap has been taken into account, the theory can be 
reduced to an effective 2x2 Hamiltonian for each sur- 
face, each of which contains a a z term due to interlayer 
tunneling. This mass term turns the original Dirac cones 
of the surface states into two massive parabolas. Twofold 
Kramers degeneracy is preserved since time-reversal sym- 
metry is not broken. Consequently, the parabolas are in- 
dexed by t z = ±1, describing two sets of Dirac fermions 
(rather than two independent surfaces). As a result of 
the mass term the Berry curvature of the surface states 
is finite, as is the orbital magnetic moment arising from 
the self-rotation of Bloch electron wave packets. The 
mass term oscillates as a function of film thickness and 
changes sign at critical thicknesses, where a topological 
phase transition occurs due to a discontinuous change in 
Chern number, in the same fashion as in the 2D quantum 
spin-Hall insulator. 

Li et studied transport in topological insulators 
in the neighborhood of a helical spin density wave. The 
nontrivial magnetic structure results in a spatially local- 
ized breakdown of time reversal symmetry on the TI sur- 
face. The helical spin density wave has the same effect as 
an external one-dimensional periodic magnetic potential 
(of period L) for spins on the surface of the TI, with an 
additional term in the Hamiltonian 



U (x) — UyCJy cos 



27T2; 



U z a z sin I — r- I . 



(23) 



As expected the Dirac cone of the TI surfaces becomes 
highly anisotropic, and, as a result, transport due to the 
topological surface states displays a strong anisotropy. 
The helical spin density wave on the TI surface leads to 
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FIG. 4: Spatial distribution of local spin density S z for a 
single Dirac fermion with a Fermi velocity of 5 x 10 5 ms -1 . (a) 
n— n junction, V2 = 40 meV, and the incoming electron energy 
is 60 meV. The angle of incidence 9 = 63 deg [(blue) solid line] 
and 8 — 18 deg [(red) dashed line], (b) n—p junction, V% = 60 
meV and the incoming electron energy is 35 meV. The angle 
of incidence 6 = 63 deg [(blue) solid line] and 6 = 36 deg [(red) 
dashed line]. This figure has been adapted from Ref. 1711 



striking anisotropy of the Dirac cones and group veloc- 
ity of the surface states. The decrease in group velocity 
along the direction defined by the one-dimensional pe- 
riodic magnetic potential is twice as large as that per- 
pendicular to this direction. More importantly, the re- 
searchers found that at the Brillouin zone boundaries, the 
periodic magnetic structure gives rise to new scmi-Dirac 
points that have a linear dispersion along the direction 
defined by the magnetic potential but a quadratic disper- 
sion perpendicular to this direction. The group velocity 
of electrons at these new semi-Dirac points is predictably 
also highly anisotropic. The authors suggest that such 
a structure could be realized in a topological insulator 
grown on a multiferroic substrate. It could be the direct 
manifestation of the elusive 2D surface TI transport even 
in the presence of considerable bulk conduction, since the 
bulk transport is isotropic and presumably unaffected by 
the presence of the helical spin density wave. 

A fascinating publication recently predicted a giant 
polarization due to the topological surface states of 
Bi2Se3p3 The structure proposed by the authors envis- 
ages a TI thin film to which a step-function potential 
is applied by means of a back gate, given by V(x) = 
V 2 Q(x). A voltage V y across the surface drives a charge 
current. The charge current in turn generates a spin- 
Hall current along the ^-direction, which strikes the step- 
potential boundary at x = 0. Near the potential bound- 
ary at x = 0, a spin-z polarization (20%) exists along the 
^-direction, induced by an electric current along the y- 



direction in the ballistic (disorder-free) transport regime 
at zero temperature. Thanks to the Klein paradox, elec- 
trons with energies less than V 2 are not confined. How- 
ever, the incident electrons are no longer in eigenstates 
on the x > side of the potential boundary. The spin po- 
larization oscillates across the potential boundary with a 
period given by 1/kp, reflecting the interference of elec- 
tron and hole states. This is shown in Fig. |4j Note that 
the plot represents the spin polarization of a single elec- 
tron. The local spin density has a strong dependence on 
the angle of incidence of the electrons. The induced spin 
polarization is found to be insensitive to the surface e^, 
because the spin polarization is approximately inversely 
proportional to the Fermi velocity, which is a constant 
for Dirac electrons. The authors expect the effect to be 
observable in a spin resolved STM. 

Still on the topic of spin transport, Dora and Moess- 
ner studied the spin-Hall effect of massive Dirac fermions 
in topological insulators in strong electric fields, near the 
threshold of electrical breakdown P23 It is well known that 
a single Dirac cone with a mass term gives rise to a quan- 
tized spin-Hall effect, as discussed above in the case of 
Ref. |69l Dora and Moessner consider switching on an 
electric field, following which a steady-state Hall current 
is reached via damped oscillations of the transient com- 
ponents. They find that the spin-Hall conductivity re- 
mains quantized as long as the electric field does not in- 
duce Landau-Zener transitions. However, as the Landau- 
Zener threshold is approached, the spin-Hall conduc- 
tance quantization breaks down, and the conductivity 
decreases asl/vE. In other words, electrical breakdown 
affects not only longitudinal transport, but also trans- 
verse transport. 

The topological order present in TI is a result of one- 
particle physics. Interaction effects have been stud- 
ied in TI with a focus on phenomena in equilibrium 
and in the quantum Hall regime. A fundamental ques- 
tion is whether basic TI phenomenology survives in- 
teractions out of equilibrium (i.e. in an electric field), 
since in transport topology only protects against back- 
scattering. A second fundamental question concerns the 
non-equilibrium spin polarization due to spin-momentum 
locking in an electric field. In interacting TI, the mean- 
field Hamiltonian is spin-dependent, and it is natural 
to ask whether a Stoner criterion exists for the non- 
equilibrium spin polarization. A systematic understand- 
ing of electron-electron interactions in non- equilibrium 
TI remains to be constructed. It was demonstrated 
recently^ that many-body interactions enhance the 
charge conductivity and non-equilibrium spin polariza- 
tion of TI, and that the non-equilibrium enhancement 
is intimately linked to Zitterbewegung and not diver- 
gent under any circumstances. Although disorder renor- 
malizations give a nontrivial doubling of the enhance- 
ment, the equivalent of the Stoner criterion is never 
fulfilled for TI. These findings can be understood by 
means of the following physical picture, illustrated in 
Fig. [5j The enhancement of the charge conductivity and 



12 




No E E // x 
> 

FIG. 5: Spin polarization in interacting TI. At each fc, red 
arrows represent the effective spin-orbit field Clk and green 
arrows the spin, (a) Equilibrium non-interacting case. The 
spin follows the local fife and the spin polarization averages 
to zero. Out of equilibrium (not shown) the spin at fc also 
follows flk, resulting in a net spin polarization, (b) Non- 
equilibrium interacting case. Interactions tilt the spin away 
from fife (here the tilt || y) and enhance the non-equilibrium 
spin polarization. Scattering from fc = fci to k — kx further 
increases the local (i.e. fc-dependent) spin polarization, while 
scattering from fc = fci to fc = — kx decreases it. In an electric 
field the asymmetry of the Fermi surface means that the net 
effect of scattering is to increase the spin polarization. This 
figure is adapted from Ref. 1731 



non-equilibrium spin polarization reflects the interplay of 
spin-momentum locking and many-body correlations. A 
spin at fc feels the effect of two competing interactions. 
The Coulomb interaction between Bloch electrons with 
k and fc' tends to align a spin at k with the spin at fc', 
equivalent to a i-rotation. The total mean-field inter- 
action tends to align the spin at fc with the sum of all 
spins at all fc', and the effective (Hartree-Fock) Hamil- 
tonian describing electron-electron interactions encapsu- 
lates the amount by which the spin at fc is tilted as a 
result of the mean-field interaction with all other spins 
on the Fermi surface. The effective field fi^ tends to 
align the spin with itself. As a result of this latter fact, 
out of equilibrium, an electrically-induced spin polariza- 
tion is already found in the non-interacting systemp2l as 
discussed extensively above. Interactions tend to align 
electron spins with the existing polarization. The effec- 
tive z-rotation explains the counterintuitive observation 
that the enhancement is related to Zitterbewegung, orig- 
inating as it does in Sk±- Many-body interactions give 
an effective fc-dependent magnetic field || z, such that for 
E || x the spins s y and — s v are rotated in opposite di- 
rections, reinforcing each other. Due to spin-momentum 
locking, a tilt in the spin becomes a tilt in the wave vec- 
tor, increasing the conductivity: spin dynamics create a 
feedback effect on charge transport. This feedback effect 



is even clearer in the fact that scattering doubles the en- 
hancement, a feature that was demonstrated rigorously 
using the density matrix formalism. This doubling is 
valid for any elastic scattering. Take fc = fc\ in Fig. [H] 
Scattering from fc = fc x to fc = kx increases the spin 
polarization at that point on the Fermi surface, while 
scattering from fc = fci to fc = —kx decreases the spin 
polarization at that point. Because the Fermi surface is 
shifted by the electric field, the net effect of this mecha- 
nism is to increase the spin polarization. 

This concludes the discussion of the theory of non- 
magnetic surface transport in topological insulators. 
Aside from the predictions of some fascinating effects, 
the central message of this section is that three primary 
signatures of surface transport may be identified: (a) 
a non-universal minimum conductivity as in graphene, 
due to the formation of electron and hole puddles as the 
Dirac point is approached (b) a current-induced spin po- 
larization and (c) a thickness independent conductance 
would be a likely indication of TI surface transport. Al- 
though not elaborated upon, the latter fact is obvious, 
since the number of surface carriers does not change with 
film thickness, wheres the number of bulk carriers does. 

IV. MAGNETOTRANSPORT THEORY 

Both in theory and in experiment, considerably more 
attention has been devoted to magnetotransport than to 
non-magnetic transport. The primary reason for this fo- 
cus is that specific signatures of surface transport can 
be extracted without having to tune through the Dirac 
point or measure the spin polarization of a current. In 
Shubnikov-de Haas oscillations and in the anomalous 
Hall effect, both discussed below, qualitatively different 
behavior is expected from the bulk and surface states. 

Broadly speaking, magnetotransport encompasses two 
categories: transport in magnetic fields and transport in 
the presence of a magnetization. A magnetic field couples 
both to the orbital degree of freedom, giving rise to the 
Lorentz force and cyclotron orbits, and to the spin degree 
of freedom via the Zeeman interaction. A magnetization 
couples only to the spin degree of freedom. The electrical 
response in the presence of a magnetic field and that in 
the presence of a magnetization do share common terms, 
yet they ultimately contain different contributions and 
are best treated separately. I will concentrate first on 
transport in a magnetic field, then on transport in the 
presence of a magnetization. Well-known aspects that 
have been covered extensively before, primarily the topo- 
logical magnetoelectric effect will not be covered here. 

A. Magnetic fields 

In a magnetic field B the energy spectra of both the 
bulk and the surface consist of Landau levels. The cy- 
clotron frequency lj c = eB/m*, with m* the effective 
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mass, is conventionally used to differentiate weak and 
strong magnetic fields. In a weak magnetic field uj c t <C 1, 
while in a strong magnetic field lo c t 3> 1, with r the mo- 
mentum relaxation time as in the previous section. In 
the bulk case, with B || z, the diagonal conductivity a xx 
contains a classical contribution as well as an oscillatory 
term periodic in 1/B, referred to as Shubnikov-de Haas 
(SdH) oscillations. Unlike non-magnetic transport, the 
classical term represents scattering- assisted transport,^ 
which can be interpreted as the hopping of electrons be- 
tween the centers of the classical cyclotron orbits. In 
both weak and strong magnetic fields this contribution 
is oc 1/t. For details of its form in different regimes 
the reader is referred to the comprehensive discussion in 
Ref. [5H1 The oscillatory contribution to the diagonal con- 
ductivity reflects the fact that consecutive Landau levels 
pass through the Fermi energy as the magnetic field is 
swept. The oscillatory contribution is present even at 
small magnetic fields, that is, in the ordinary Hall effect. 
The corresponding oscillations in the magnetization are 
termed the de Haas- van Alphen (dHvA) effect. In the 
bulk in a magnetic field only the ordinary Hall effect will 
occur, and experimentally one expects SdH oscillations 
for all orientations of the magnetic field. 

In a strong perpendicular magnetic field in 2D systems 
the quantum Hall effect is observed. The diagonal con- 
ductivity a xx consists of a series of peaks, while the Hall 
conductivity a xy consists of a number of plateux. The 
Landau levels corresponding to the topological surface 
states are effectively the same as in graphene and the 
theory will not be developed here, primarily be caus e it 
is either obvious or it has been covered elsewhere Ne- 
glecting the Zeeman splitting, the Landau level energies 
indexed by n are given by the formula 

e n = e + sgn(n) - v / 2heB\nl (24) 

with Sq the energy at the Dirac point, which is inde- 
pendent of the magnetic field for massl ess R ashba-Dirac 
fermions. The Landau levels show the yJ\n\B dependence 
on the magnetic field and Landau level index character- 
istic of the linear Rashba-Dirac cone dispersion. In two- 
dimensional systems the considerations above for a xx re- 
main valid, while in a large perpendicular magnetic field 
the quantization of a xy becomes apparent. The quan- 
tized Hall conductivity is given by 

<T xy = (n+ G —. (25) 

The half-quantization is due to the 7r Berry phase, which 
is also present in the scattering term in non-magnetic 
transport. The half-quantization of the QHE can also be 
understood in terms of a term in the bulk Lagrangian 
oc E ■ B. Shubnikov-de Haas oscillations in a xx are de- 
scribed by the same formalism as in graphene. Only the 
levels at the Fermi energy at one particular time con- 
tribute to the diagonal conductivity cr xx , yet the Hall 



conductivity a xy has contributions from all levels below 
the Fermi energy. 

In a magnetic field of several Tesla, if the magneto- 
transport signal can be ascribed to the topological sur- 
face states, one thus expects to observe the quantum Hall 
effect, and expects to see SdH oscillations only for a mag- 
netic field normal to the TI surface. An in-plane mag- 
netic field would not yield SdH oscillations. In addition, 
due to the different effective mass, the frequency of SdH 
oscillations for the bulk is considerably different from 
that for the surface. Specifically, the cyclotron frequency 
u> c = eB/m* is a constant for bulk transport. For sur- 
face transport one may define a cyclotron effective mass 
m* = h 2 kp/A, which depends on the Fermi wave vector, 
given by kp — \fAim. In the case of coexisting bulk and 
surface transport one therefore expects two sets of SdH 
oscillations with different frequencies, as well as angular 
dependence of the SdH oscillations, both of which should 
be clearly distinguishable experimentally. 

Shubnikov-de Haas and de Haas-van Alphen oscilla- 
tions in Bi 2 Se3 were studied theoretically in Ref. 1741 The 
authors determine the band structure from first princi- 
ples, then calculate the magnetization induced by an ex- 
ternal magnetic field. Depending on whether the zero- 
mode Landau level is occupied or empty, the location 
of center of the dHvA oscillations in the magnetization 
changes. Interestingly, the magnetic oscillation patterns 
for the cases of filled and half-filled n — Landau level 
are out phase. 

Chu et aZP^ studied numerically the transport proper- 
ties of the topological surface states in a Zeeman field. 
The main argument of their work is that a chiral edge 
state forms on the surface, but is split into two spatially- 
separated halves. Each half carries one half of the con- 
ductance quantum. They propose a four-terminal setup 
in which this theory can be tested. Numerical simula- 
tions reveal that the difference between the clockwise and 
counterclockwise transmission coefficients of two neigh- 
boring terminals is ~ 1/2. This suggests a half-quantized 
ordinary Hall conductance could be measured in a four- 
terminal experiment. 

Zhang et aZP^set out to explain why the half-quantized 
QHE has not been observed in experiments to date. The 
authors investigated a 3DTI in a uniform magnetic field 
using tight binding model. They point out that all surface 
states of a finite 3DTI are connected to each other. In 
a finite 3D TI, the two surfaces cannot simply be sepa- 
rated by the bulk as in usual 2D quantum Hall systems. 
As a result, surface states may move from one side of 
the sample to another. Hall voltage measurements for 
3DTI are therefore necessarily ambiguous. The Hall cur- 
rent is due to 2D surface states, rather than ID edge 
states. The structures of these side surface states are 
more complicated than the edge states in 2DEG. Quan- 
tum Hall plateaux of (2n + 1) e 2 /h can only be defined 
from transverse current instead of the Hall voltage. The 
tight-binding calculations are in agreement with earlier 
predictions based on an effective theory of Dirac fermion 
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in curved 2D spaces,^ which showed that, based on the 
argument above that the current is carried by 2D rather 
than ID states, the robustness of the QH conductance 
against impurity scattering is determined by the oddness 
and evenness of the Dirac cone number. 

Tse and MacDonalcP^ developed a theory of the 
magneto-optical and magneto-electric response of a topo- 
logical insulator thin film in a strong (quantizing) mag- 
netic field. Their work was driven by the realization that 
the quantum Hall effect due to the topological surface 
states could be observed optically, even when bulk carri- 
ers are present in the system and contribute to conduc- 
tion. Interestingly, these authors showed that the low- 
frequency magneto-optical properties depend only on the 
sum of the Dirac-cone filling factors of the top and bot- 
tom surfaces. In contrast, the low-frequency magneto- 
electric response is determined by the difference in these 
filling factors. Furthermore, the Faraday rotation angle 
is quantized, being given by the sum of the filling factors 
of the top and bottom surfaces multiplied by the fine 
structure constant. The Faraday rotation angle exhibits 
sharp cyclotron resonance peaks and changes sign near 
each allowed transition frequency, as shown in Fig. [6] 
At half-odd-integer filling factors a single dipole-allowed 
intraband transition exists. At other filling factors two 
resonances are allows at different frequencies associated 
with transitions into and out of the partially filled Lan- 
dau level. This behavior is completely different from 
usual 2DEGs, in which all dipole-allowed transitions have 
the same energy. The Kerr rotation angle is 7r/2, as 
found by the same authors in the case of magnetically- 
doped TI (see below.) It is strongly enhanced at the 
frequencies at which cyclotron resonance occurs. Unlike 
the low- frequency giant Kerr effect, the resonant peaks 
corresponding to absorptive Landau level transitions are 
non-universal. 

Tkachov and Hankiewicz^ studied the quantum Hall 
effect due to the topological surface states in a TI film 
of thickness d in the situation in which bulk carriers also 
contribute via the ordinary Hall effect. The authors took 
into account the contributions of the bulk as well as the 
top and bottom surfaces, showing that the contribution 
to the conductivity due the bulk carriers can be sup- 
pressed by an external magnetic field and an ac-electric 
field. They work in the regime uj,lo c <C A, where u> is 
the frequency of the applied electric field, and A is the 
Zeeman-energy gap in the Dirac cone. Under these con- 
ditions, the cyclotron resonance frequency is shifted by a 
thickness-dependent term, 



da' 



bulk 
xx 
sfc\ 



(26) 



where o^x is the bulk dc conductivity at zero mag- 
netic field, and a s Jy is the surface Hall conductivity. The 
nonzero surface conductivity alters the cyclotron drift in 
the direction of the current, and induces a linear nega- 
tive magnetoresistivity. At the same time, the Hall angle 
depends on the magnetic field quadratically 
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FIG. 6: (a) Faraday rotation angle versus frequency oj/(ec) 
at equal filling factors on both surfaces for a 30nm-thick 
Bi2Se3 film. The densities on both surfaces are 5 x 10 n cm~ 2 . 
The bulk band gap is 0.35 eV, the Fermi velocity is up = 
5 x 10 ms _1 and the dielectric constant e r — 29. (b). Kerr 
rotation angle versus frequency for the same parameters. This 
figure is adapted from Ref. [78] 



Zyuzin et alW^ investigated a thin TI film in a strong 
perpendicular magnetic field, allowing for hybridization 
between the top and bottom surfaces. The authors deter- 
mined the Landau-level spectrum of the film as a function 
of the applied magnetic field and the magnitude of the 
hybridization matrix element, taking into account both 
the orbital and the Zeeman spin splitting effects of the 
field. Most interestingly, they identified a quantum phase 
transition between a state with a zero Hall conductivity 
and a state with a quantized Hall conductivity of e /h 
as a function of the magnitude of the applied field. The 
transition is driven by the competition between the Zee- 
man and the hybridization energies. 



B. Magnetizations 

Doping with magnetic impurities, whose moments 
align beyond a certain doping concentration, makes topo- 
logical insulators magnetic and opens a gap. Topological 
insulators in the presence of a strong out-of-plane magne- 
tization exhibit rather different physics than their coun- 
terparts placed in magnetic fields. 

The interplay of magnetism, strong spin-orbit cou- 
pling, disorder scattering and driving electric fields leads 
to the anomalous Hall effect in magnetic TIs. The 
anomalous Hall effect has been a mainstay of condensed 
matter physics, studied extensively over the years.^ 
A recent groundbreaking paper by Yu et alW^ studied 
tetradymite semiconductors doped with Cr and Fe using 
first principles calculations, determining a critical tem- 
perature as high as 70K. Yu et al. further demonstrated 
that in a thin TI film near a band inversion, in which the 
chemical potential lies in the magnetization-induced gap 
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between the surface valence and conduction bands, edge 
states give rise to an integer quantized anomalous Hall 
effect The quantized anomalous Hall effect requires a 
band inversion, provided by spin-orbit coupling, and fer- 
romagnetic order, provided by doping. The 4x4 Hamil- 
tonian describing the surface states can be broken up into 
two blocks, and a sufficiently large exchange field gives 
rise to a difference in Chern number of one between the 
two blocks, resulting in an anomalous Hall conductivity 
of e 2 /h. Remarkably, this is true regardless of whether 
the system is in the topologically trivial phase or in the 
topologically insulating phase. The quantized anoma- 
lous Hall effect is to be distinguished from the anoma- 
lous quantum Hall effect, which has also been predicted 
to occur in topological insulatorsP 

In contrast to the 2D case, in 3D magnetic TIs with 
the chemical potential in the band gap the anomalous 
Hall condu ctivity due to one surface is half- quantized 
e 2 /2h} 83 ' ' rather than the integer quantization seen for 
2D edge states. The value of e 2 /2h found for Ep in 
gap was demonstrated to be robust against disorder.^ 
Ref . l85l studied the anomalous Hall effect due to the sur- 
face states of 3D topological insulators in which e p lies in 
the surface conduction band, showing that the topologi- 
cal band structure contribution (e 2 /h) found in previous 
studies, with a substantial renormalization due to dis- 
order, yields the dominant term, overwhelming all con- 
tributions due to skew scattering and side jump. The 
end result is of the order of the conductivity quantum 
and independent of the magnetization. These results are 
explained in what follows. 

Theoretically, the already challenging problem of the 
anomalous Hall effect in 3D TI is complicated by the ne- 
cessity of a quantum mechanical starting point and of 
a matrix formulation describing strong band structure 
spin-orbit coupling, spin-dependent scattering including 
rclativistic corrections, Klein tunneling and Zitterbewe- 
gung on the same footing. The density matrix-Liouville 
equation formalism has a quantum mechanical founda- 
tion and accounts for the three mechanisms known to 
be important in the anomalous Hall effect. The inter- 
play of band structure spin-orbit interaction, adiabatic 
change in k in the external electric held, and out-of-plane 
magnetization leads to a sideways displacement of carri- 
ers even in the absence of scattering (intrinsic)W^ The 
spin dependence of the impurity potentials via exchange 
and extrinsic spin-orbit coupling causes asymmetric scat- 
tering of up and down spins (skew scattering),^ and a 
sideways displacement during scattering (side jump)P^ 

In magnetically doped TI the magnetization M ^ 
and the kinetic equation can be written in the forr 



,85 



dfk 
dt 



[H m ki fk] + J(fk) = T^k-, 



with H m k given by Eq. Q . All extrinsic spin-orbit terms 
oc A must be taken into account, including the sec ond 
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Born approximation scattering term J (fk) of Eq. 
the spin-dependent correction to the position operator 



oc A and the corresponding modification of the interac- 
tion with the electric field. Consequently, the driving 
term T>f. contains a host of contributions not found in 
the non-magnetic case. It is given in detail in Ref. 1851 
At zero temperature the splitting due to M is resolved, 
and the Fermi energy > \M\ as usual lies in the sur- 
face conduction band. Three orthogonal directions are 
identified in reciprocal space, 



&k = -ak + b k z 



k efJ = k 



(28) 



J e// — a k Z + Ok 



he, 



where a k = 2Ak/Mlk and b k — 2M/Mlk, with fi^ de- 
fined in Eq. Q, and ap and &f at k — hp, with bp -C 1 
at usual transport densities. Considering the negligible 
size of 6fc at k — kp, the condition yielding the suppres- 
sion of backscattering is effectively unmodified from the 
non-magnetic case. The vector cr of Pauli matrices is 
projected onto the three directions in reciprocal space in 
an analogous fashion to the non-magnetic case, and the 
details will not be reproduced here .^Likewise Sk is also 
projected onto these three directions. 

The kinetic equation is solved perturbatively in 
h/(epr), characterizing the disorder strength, as well as 
in A, which quantifies the strength of the extrinsic spin- 
orbit coupling. When a magnetization is present a num- 
ber of terms independent of r appear in the conductivity, 
representing transvers (Hall) transport. The dominant 
term in the anomalous Hall conductivity is^ 



2h 



(1-a), 



(29) 



where a is a disorder renormalization, equivalent to a 
vertex correction in the diagrammatic formalism.^ The 



bare contribution to a yx is 



2h' 



which can also be ex- 



pressed in terms of the Berry (geometrical) curvature^, 
and is thus a topological quantity. It is similar to the re- 
sult of Ref. 83 in the vicinity of a ferromagnetic layer, and 
can be identified with a monopolc in reciprocal space. It 
arises from the integral 



AkM 



(A 2 k 2 + M 2 ) 3 / 2 



M 



y/A 2 k' 2 p + M 2 



(30) 



If we ignore terms of order b 2 F it can be approximated 
as simply 1/A, corresponding to the limit kp — > oo. In 
this regime the integral is therefore a constant, and the 



(27) bare conductivity a yx 



2h 



independent of the magne- 



tization, number density, or Rashba spin-orbit coupling 
strength A. Contributions oc A turn out to be negligible 
in comparison with this term.^ 

Ref. 85 showed that there are no terms from extrin- 
sic spin-orbit scattering that commute with the band 
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Hamiltonian, and thus, following the argument above, 
no terms in the density matrix oc r due to extrinsic 
spin-orbit mechanisms. Neither extrinsic spin-orbit cou- 
pling nor magnetic impurity scattering give a driving 
term in the kinetic equation parallel to Hq^ (that is, 
|| fife), thus the anomalous Hall response does not con- 
tain a term of order [h/ {eft)]~ 1 due to extrinsic spin- 
orbit scattering or magnetic impurity scattering. Such a 
term would have overwhelmed the renormalized topolog- 
e 2 

ical term <t„ t = (1 — a) in the ballistic limit. 

In TI the spin and charge degrees of freedom are in- 
herently coupled and the anomalous Hall current can 
also be viewed a steady-state in-plane spin polarization 
in a direction parallel to the electric field. The (bare) 
topological term has two equal contributions, which are 
part of the correction to the density matrix orthogonal 
to the effective Zeeman field, therefore they represent an 
electric-field induced displacement of the spin in a direc- 
tion transverse to its original direction. They are also ob- 
tained in the Heisenberg equation of motion if one takes 
into account the fact that k is changing adiabatically. 
Physically, the ^-component of the effective Zeeman field 
Jlfc is changing adiabatically, and the out-of-plane spin 
component undergoes a small rotation about this new ef- 
fective field. Consequently, each spin acquires a steady 
state component parallel to E, which in turn causes k to 
acquire a small component in the direction perpendicular 
to E. Elastic, pure momentum scattering [contained in 
Joifk)] reduces this spin polarization because, in scatter- 
ing from one point on the Fermi surface to another, the 
spin has to line up with a different effective field Jl^. The 
extra spin component of each electron is tx M, however 
the final result is independent of M, as the integr and co n- 
tains a monopole located at the origin in fc-spacel- 83 * 89 * As 
M — > the effect disappears, since the correction to the 
orthogonal part of the density matrix vanishes. Finally, 
though the form of this term would hint that it is observ- 
able for infinitesimally small M, it was tacitly assumed 
that M exceeds the disorder and thermal broadening. 

The appearance of terms to order zero in h/(eFT), 
not related to weak localization, is standard in spin- 
dependent transport. The expansion of the density ma- 
trix starts at order h/(eFT'Y~ 1 \ i.e. the leading-order 
term is oc r, so the fact that the next term is indepen- 
dent of /i/(ei?r)(~ 1 ) is to be expected. Of course this 
term is not independent of the angular characteristics of 
the scattering potential (a above depends on the scatter- 
ing potential), which is why it is not enough to consider 
only the band structure contribution, but also its disor- 
der renormalization. Disorder renormalizations are cru- 
cial in spin-related transport, where they can go as far as 
to cancel band structure contributions, as they do in the 
spin-Hall effect in spin-orbit coupled systems described 
by the Rashba Hamiltonian P3 

Experimentally, for charged impurity scattering the 
figures depend on the Wigner-Seitz radius r si represent- 
ing the ratio of the Coulomb interaction energy and the 



kinetic energy, with e r the relative permittivity. For 
Bi2Se3 with r s = O.l^P^l, the dominant term by far is 
o~yx ~ — 0.53 (e 2 /2h) « — e 2 /4h. We can consider also 
the (artificial) limit r s — > 0, which implies e r — > oo, arti- 
ficially turning off the Coulomb interaction. As r s — > 0, 
the prefactor of —e 2 /2h tends to 0.61, while at r s = 4, 
the limit of RPA in this case, it is ~ 0.12. Interestingly, 
for short-range scattering, the anomalous Hall current 
changes sign, with a yx ~ 0.18 (e 2 /2h). The expression 
above represents the contribution from the conduction 
band, and in principle an extra e 2 /2h needs to be added 
to the total result to obtain the signal expected in ex- 
periment. Since the conduction band is expected to con- 
tribute sa —e 2 /Ah, this does not make a difference in ab- 
solute terms. However, if one bears in mind that the bare 
contribution from the conduction band is —e 2 /2h, the 
inescapable conclusion is that the remaining, observable 
part of the conductivity, is exactly the disorder renormal- 
ization a. The same issue as in the quantum Hall effect 
is present with regard to the detection of the anomalous 
Hall effect, since there is always more than one surface 
and the surfaces are connected. Unlike the quantum Hall 
effect, the surfaces give opposite contributions to the half- 
quantized anomalous Hall effect ™ But since the disorder 
renormalizations should be different on the two surfaces 
a finite total signal should exist. 

In a related context, Tse and MacDonalcP^ studied a 
TI thin film weakly coupled to a ferromagnet, but fo- 
cused on the magneto-optical Kerr and Faraday effects, 
since these are less likely to be affected by unintentional 
bulk carriers. Using linear-response theory, these scien- 
tists discovered that at low frequencies the Faraday rota- 
tion angle has a universal value, which is determined by 
the vacuum fine structure constant, when ef lies in the 
Dirac gap for both surfaces. The Kerr rotation angle in 
this regime, is again universal and equal to 7r/2, hence 
the nomenclature of a giant Kerr rotation. The physics 
overlaps significantly with the anomalous Hall effect de- 
scribed above leading to the half-quantized anomalous 
Hall effect, except in this case the effect explicitly involves 
both surfaces. The effect reflects the interplay between 
the chiral nature of the topological surface states and 
interference between waves reflected off the top surface 
and waves scattered off the bottom surface. The surface 
Hall response creates a splitting between the reflected 
left-handed and right-handed circularly polarized fields 
along the transverse direction. The presence of the bot- 
tom topological surface is important since the first-order 
partial wave reflected from the top surface undergoes only 
a small rotation. The higher-order partial waves from 
subsequent scattering with the bottom surface yield a 
contribution that strongly suppresses the first-order par- 
tial wave along the incident polarization plane. As a 
result, the left-handed and right-handed circularly polar- 
ized fields each acquire a phase of approximately 7r/2 in 
opposite directions, leading to a 7r/2 Kerr rotation. 

This concludes the discussion of the anomalous Hall 
effect. The topic of spin transfer in magnetic topological 
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insulators has also been studied recently. Nomura and 
Nagaosa^l investigated magnetic textures in a ferromag- 
netic thin film deposited on a 3D TI, showing that they 
can be electrically charged thanks to the proximity effect 
with the Dirac surface states. They consider an insulat- 
ing magnet on the TI surface, with a magnetization M 
pointing along. Due to the interplay of the chiral surface 
states with the proximity-induced magnetization, electric 
charge and current densities are induced in the TI, 



Ma 



H 



V • M 



evp 
Ma H \ dM 



evF 



(31) 



dt 



where an is the quantized Hall conductivity. The electri- 
cal current, which we have established corresponds to a 
steady-state spin density, is the analog of the spin trans- 
fer torques. Spin transfer in systems with strong spin- 
orbit interactions was discussed in the related context of 
ferromagnetic semiconductors in Ref. 93. In topological 
insulators, precisely because of the equivalence of charge 
currents and spin densities, this electric charging of mag- 
netic textures can be used for electrical manipulation of 
domain walls and vortices. In the case of vortices, an 
electric field breaks the vortcx-antivortex symmetry on 
the TI surface. In the case of domain walls, the authors 
show that an electric field can displace a domain wall in 
a topological insulator, much like a magnetic field can 
accomplish this in an ordinary magnet. 

The action of an electrical current on a magnetization 
at the interface between a topological insulator and a fer- 
romagnet was also studied by Garate and Franz P^I While 
their results were similar, these authors focused on the 
magnetization dynamics induced by the steady-state spin 
polarization, inlcuding the possibility of current-induced 
magnetization reversal. In their language, the Hall cur- 
rent changes the effective anisotropy field, which in turn 
affects the magnetization dynamics in the ferromagnet. 
Interestingly, their work shows that the magnetization 
may be flipped by the Hall charge current in tandem 
with a small magnetic field of approximately 20mT. 

The same authors subsequently focused on the magne- 
toelectric response of a disorder-free time-reversal invari- 
ant TlP^ By determining the spin-charge response func- 
tion, these researchers confirmed the observation that 
a uniform static in-plane magnetic field does not pro- 
duce a spin polarization in topological insulators. The 
equilibrium RKKY interaction between localized spins 
separated by R decays as 1/i? 3 for kp — 0, while for 
kpR ^ 1 it oscillates and decays as 1/i? 2 . Unlike ordi- 
nary two-dimensional electron systems, the RKKY coef- 
ficients depend on the carrier number density. The equi- 
librium RKKY interaction can be strongly modified by 
moderate electric fields, in particular when the localized 
impurities are separated by distances in excess of 5nm. 
The authors also found that the helical metal mediates 
a strong Dzyaloshinskii-Moriya coupling at finite doping, 



which favors noncollinear magnetization configurations. 
Lastly, thanks to the inherent spin-momentum locking 
of topological insulators, an impurity magnetic moment 
that precesses in an external Zeeman field gives rise to 
an alternating charge current. 

Other researchers have brought to light different facets 
of the exciting new physics of magnetic topological insu- 
lator. A recent work by Burkov and Hawthorn 96 focused 
on spin-charge coupling due to the topological surface 
states of a diffusive topological insulator. The authors 
derived a set of diffusion equations describing spin-charge 
coupled transport of the helical metal surface. The most 
interesting finding of Ref. [96] is that the continuity equa- 
tion for the carrier density n and particle current density 
j a new, spin-dependent term appears, 
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dt 
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-DVn — vp S x z, 



where D is the usual diffusion coefficient and S represents 
the spin density. The second, unconventional term in the 
current, related to the chirality (helicity) of the topo- 
logical surface states, leads to a new magnetoresistance 
effect. The authors envisage a TI surface on which cur- 
rent flows between a ferromagnetic electrode and a non- 
magnetic electrode. The ferromagnetic electrode injects 
a spin-polarized current. Solving the diffusion equation, 
it emerges that the voltage drop consists of two contri- 
butions, the first being the usual Ohmic resistance, while 
the second is given by 
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(33) 



with I the current and the dimensionless rj parameteriz- 
ing the spin polarization of the current. This contribution 
is proportional to the cross product of the polarization 
of the ferromagnetic electrode and the charge current, 
and disappears when the spin polarization of the ferro- 
magnetic electrode is along the transport direction. The 
effect is independent of disorder and of the separation be- 
tween the electrodes and only depends on the spin polar- 
ization of the injected current and the Fermi momentum. 
However, the degree of the injected spin polarization de- 
pends on the (non-universal) properties of the interface 
between the ferromagnetic electrode and the topological 
insulator. The authors further show that this effect is 
tunable by a gate voltage, enabling operation of the de- 
vice as a transistor. 

Ref. [97] analyzed the tunneling conductance of a series 
of ballistic structures involving topological insulators and 
quantizing magnetic fields or magnetic materials. For a 
TI in crossed electric and strong (quantizing) magnetic 
fields, with the magnetic field perpendicular to the sur- 
face, the orbital and Zeeman effects due to an applied 
magnetic field can be modulated by the electric field. 
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The Landau levels have the form 

j? ru ^ j- hvF /oi i , 9 2 ^ 2 bJB sm6 

E ^ k y) = ± J B ^j2 W 2 M + P hv F k v> 

(for n^O) 

\9VbB sin 9\ — phvpk yi (for n = 0). 

7 

(34) 

Here the magnetic length = \/^/ (e_B sin.0), the boost 
f3 = E/(v F BsmO) < 1, and 7 = (1 - p 2 )' 1 / 2 . This 
equation shows that the Landau levels can be collapsed 
by varying either the strength of the electric field E or 
the out-of-plane tilt 9 of the magnetic field. The tunnel- 
ing conductance depends on 9, a fact easily understood 
by recalling that the in-plane component of the magnetic 
field can be removed from the Dirac-cone by a gauge 
transformation, and is therefore irrelevant. The Zeeman 
term yields an exponential suppression of the tunneling 
conductance, which is explained by the fact that the Zee- 
man term in TI is equivalent to a mass term. In nor- 
mal metal-magnetic film-normal metal (NMN) junctions 
transport can be controlled using the proximity-induced 
exchange field due to a ferromagnetic film. At a critical 
value of this exchange field, the tunneling conductance of 
the junction changes from oscillatory to a monotonically 
decreasing function of the junction width d. In principle 
such a setup could realize a magnetic switch. Normal- 
metal-barrier-magnetic film (NBM) junctions are char- 
acterized by the height Vq of the barrier. In the ab- 
sence of a magnetization, such a junction displays max- 
ima in the tunneling conductance at integer quantized 
values n of the quantity eVod/hvp. In the magnetic case, 
beyond a critical value of the exchange field, trasnmis- 
sion resonances occur at values n + 1/2 of this quan- 
tity. At even higher exchange fields the junction con- 
ductance can be fully suppressed. For normal-metal- 
magnetic film-superconductor (NMS) junctions, the posi- 
tion of the peaks of the zero-bias tunneling conductance 
can be tuned using the magnetization of the ferromag- 
netic film. 

Yokoyama et alW^ investigated conduction across a fer- 
romagnet /feromagnet interface on the surface of a topo- 
logical insulator. The device considered is the analog of 
the more usual spin-valve devices, except the substrate 
in this case is a topological insulator (the current itself 
flows in the TI) . The tunneling conductance across the in- 
terface has a strong dependence on the directions of the 
magnetizations of the two ferromagnets, related to the 
manner in which the wavefunctions on the two sides of 
the interface connect. A misfit of the in-plane momentum 
between the two sides gives rise to a strong dependence of 
the tunneling conductance on the in-plane rotation angle 
<p. Unlike conventional spin-valve devices, the conduc- 
tance in the antiparallel configuration can greatly exceed 
that in the parallel configuration. A representative part 
of their results is reproduced in Fig. [7J 

Magnetic impurity scattering has been considered in 
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FIG. 7: Anomalous magnetoresistance for a ferrromagnet- 
ferromagnet interface on top of a topological insulator. 
The figure shows the tunneling conductance for the case 
when the barrier potential between the two regions is 
zero. On the incident side the magnetization is defined as 
mi (sin 9 cos (p, sin 9 sin (j>, cos 9) while on the other side of the 
interface it is (0, 0, m^). In (a) and (c) below — 0, while in 
(b) and (d) n%2 is finite. In (a) and (b) the results are plotted 
for an n — n junction, while in (c) and (d) they are plotted 
for a p — n junction. The figure is adapted from Ref. 1981 



TI transport j 84 | 8S | although work to date solely focused 
on straightforward potential scattering, without consid- 
ering in detail spin-flip terms that lead to the Kondo 
effect. All studies above require the magnetization to be 
large enough so that it is resolved, implying that the local 
moments are all lined up, thus neglecting Kondo physics 
in the anomalous Hall effect seems sensible. Transport 
theory can be modified to include spin-flip scattering and 
Kondo physics, provided that scattering terms cx a x ,a y 
are allowed, and provided that one works explicitly in the 
grand canonical ensemble in order to account fully for 
the spin-flip factors. In the context of the Kondo effect, 
local moment formation and Friedel oscillations associ- 
ated with magnetic impurities in TI have been studied in 
Ref. [SH In addition, Biswas and Balatsky2221 calculated 
the t-matrix of TI, determining the modication of the lo- 
cal electronic structure caused by a single local impurity 
on the surface of a 3DTI. They found that the local den- 
sity of states around the Dirac point is disrupted near the 
impurity by the creation of low-energy resonance states, 
yet this does not destroy the Dirac point locally. The 
RKKY coupling between spin textures created near the 
magnetic impurities was found to be anisotropic. Aside 
from these papers, the study of the Kondo effect in TI 
is in its infancy, and there is no reliable indication of 
how high the Kondo temperature is expected to be, and 
the nontrivial effect of spin-momentum locking on Kondo 
singlet formation has not been determined. 
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V. WEAK LOCALIZATION 

On a fundamental level, weak localization is an inter- 
ference phenomenon between two time-reversed paths in 
a scattering event, which gives a correction to the con- 
ductivity of a solid. Weak localization refers to interfer- 
ence between the incident and backscattered paths, and 
is often referred to as coherent backscattering. To obtain 
the weak localization correction in the linear response 
Kubo formalism, one would have to return to Eq. [9] and 
determine the impurity average of the product of two 
Green's functions in the next order in the small parame- 
ter H/(eFT). Logistically, this is tantamount to calculat- 
ing the Cooperon in the Kubo formula, that is, the sum 
of maximally crossed diagrams, which can be mapped 
onto a sum of twisted ladder diagrams. In the clas- 
sical interpretation of this phenomenon, the maximally 
crossed diagrams represent the probability that a parti- 
cle will return to its starting position after a scattering 
process. Although the Liouville-density matrix formal- 
ism has not been extended to cover weak localization, an 
analogous and widely used kinetic equation formalism 
can be devised starting from the Keldysh Green's func- 
tion approachPS This can be straightforwardly extended 
to incorporate weak localization, and the weak localiza- 
tion correction to the conductivity in this language is also 
expressed in terms of the Cooperon. 

Being an interference phenomenon, weak localization 
reflects the existence of phase coherence between two 
time-reversed paths, incident and backscattered. This is 
quantified by defining a phase relaxation time r^, which 
represents the time over which the electron wave function 
retains its phase. Phase relaxation is typically caused by 
dynamic perturbations such as electron-electron scatter- 
ing and phonon scattering, although in practice the phase 
relaxation time is frequently taken as a phenomenologi- 
cal quantity. In an ordinary 2DEG the weak localization 
correction to the conductivity takes the form 
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Phase coherence between time-reversed paths is also 
destroyed by magnetic interactions, thus weak localiza- 
tion is correspondingly suppressed. For example, there 
is no weak localization correction to the ordinary Hall 
effect. Magnetic impurities also reduce weak localiza- 
tion, although they do not eliminate it. The disappear- 
ance of weak localization in a magnetic field makes this 
phenomenon a useful experimental characterization tool. 
Weak localization gives a negative contribution to the 
conductivity, therefore application of a magnetic field 
yields a positive magnetoconductance (negative magne- 
toresistance.) In 2D weak localization has a logarithmic 
temperature-dependence, and electron-electron interac- 
tions also lead to a conductivity correction logarithmic 
in the temperature. 

If spin-orbit coupling is present in the impurity poten- 
tial or in the band structure, the correction to the con- 



ductivity is positive, and is termed weak antilocalization. 
The same is true for the case of a pseudospin dependent 
Hamiltonian, such as that of graphene monolayers and 
bilayersp^H and is expected to apply to topological insu- 
lators. The occurrence of weak antilocalization in topo- 
logical insulators can be easily understood by recalling 
that backscattering is absent. Simply put, there cannot 
be coherent backscattering when there is no backscatter- 
ing at all. Weak antilocalization manifests itself in exper- 
iment as a negative magnetoconductance (positive mag- 
nctoresistance.) The full magnetic field dependence in 
an ordinary 2DEG is encapsulated in the Hikami-Larkin- 
NagaokcP2 formula for 5a™ 9 = 5a xx {B) - Sa xx (0) 
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where \& represents the digamma function, which is the 
logarithmic derivative of the Eulcr gamma function, and 
D = vft 2 /2 the diffusion coefficient. 

At sufficiently high impurity densities, coherent 
backscattering leads to Anderson localization, and the 
material undergoes a transition to an insulating state 
which bears the name of the Anderson metal-insulator 
transition. Although in topological insulators the ab- 
sence of backscattering is believed to suppress Ander- 
son localization, one remark is in order. This problem 
has been studied with the aid of renormalization group 
techniques, and it was determined that the lower criti- 
cal dimension for the localization problem is dim = 2. 
This seems to imply that all states are localized in 2D. 
Yet the bare textbook argument that no extended states 
exist in 2D is not borne out by experiment, since the lo- 
calization length can easily exceed the size of the sample. 
Taken literally, this would imply that the 2D electron gas, 
MOSFETs and other semiconductor devices do not exist, 
whereas quasi-2D systems display the highest known mo- 
bilities in the world (tens of millions cm 2 /Vs in GaAs), 
as a result of modulation dopi ng. 

Tkachov and HankiewiczSSSl studied weak antilocaliza- 
tion in HgTe quantum wells, which have a finite mass, 
and topological surface states, which are massless. In 
HgTe quantum wells, the weak localization correction is 
shown to be 
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where a is a constant, r e is the elastic scattering lifetime, 
while r m is given by 
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(38) 



and the mass parameter in the Dirac equation is m+Bkp, 
with m and B material specific constants, and A as de- 
fined in the Rashba-Dirac Hamiltonian. The finite mass 
suppresses weak antilocalization in HgTe quantum wells 
at times larger than r m , which corresponds to the open- 
ing of a relaxation gap h/T rn in the Cooperon diffusion 
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mode. Since in HgTe quantum wells m can be either pos- 
itive or negative depending on the well width, r m has a 
non-monotonic dependence on both m and n, and van- 
ishes in the case when 



2imB = 0. 



(39) 



Along this line the weak antilocalization term in the con- 
ductivity reaches its maximum value, given by the usual 
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In the case of topological 



insulators, the authors find that the weak antilocalization 
correction takes the form 
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the same as a 2DEG with spin-orbit interaction. How- 
ever, in a parallel magnetic field, due to the exponential 
decay of the surface states into the bulk, topological in- 
sulators exhibit a weak-antilocalization magnetoconduc- 
tivity which is a function of the penetration length into 
the bulk. The time scale for the suppression of the quan- 
tum interference by the magnetic flux is determined by 
the thickness of the surface state, which is expected to 
be much smaller than the magn etic length. 

Concomitantly, Lu et a/P^ identified a competition 
between weak localization and anti-localization in topo- 
logical surface states. They find that, as expected, the 
gapless topological surface states indeed give rise to weak 
antilocalization. At the same time, however, when the 
gap is opened (by magnetic doping, for example), a new 
weak localization contribution emerges. The gap and 
magnetic scattering can change the magnetic field de- 
pendence from weak antilocalization-like to a parabolic 
dependence, and increasing the size of the gap with re- 
spect to ep may drive the system into the weak local- 
ization regime. In this regime, however, the weak local- 
ization correction is half that found in the usual 2DEG. 
The transition from weak antilocalization to weak local- 
ization is understood as a change in Berry phase as the 
magnetization gap increases from ze ro to 2ep. 

In a subsequent paper, Lu and Shert^l concentrated on 
the bulk channel contribution to weak localization. The 
authors argue that the recently observed weak antilocal- 
ization in TI appears to have contributions from the bulk 
states. Therefore they devised a unified model that al- 
lowed them to study bulk and surface contributions to 
weak localization and antilocalization and magnetocon- 
ductivity. Within the framework of fc • p perturbation 
theory the authors derived a modified Dirac Hamiltonian 

H 3D = C+D(k 2 + k 2 z )+Aa-k+[m-B(k 2 + k 2 z )} f3, (41) 

where k — {k x , k y ), a and /3 are the usual matrices with 
the same name found in the Dirac equation, and A, B, 
C, D and m are model parameters. From this Hamil- 
tonian the topological surface bands and 2D bulk sub- 
bands are easily determined. Following a linear response 
calculation along the lines of their previous work cited 



above, the authors derived the magnetoconductivity for 
both the 2D bulk subbands and surface bands. They find 
that the bulk states contributes a weak localization term 
in the magnetoconductivity, which, in certain parame- 
ter regimes, may compensate the weak antilocalization 
terms arising from the surface states. When the bulk 
weak localization channels outnumber the surface weak 
antilocalization channels a crossover will occur at a finite 
magnetic field. The authors suggest that the competition 
between the two mechanisms may explain the observed 
experimental deviation of 5a xx from the Hikami-Larkin- 
Nagaoka formula. 

Finally, an interesting paper by Ostrovsky et aZP^ 
focused on electron-electron interactions in disordered 
topological insulators. The authors used the perturba- 
tive renormalization group to calculate the scaling func- 
tion /3(a xx ) = d1** L , where the dimensionless a xx is the 
conductivity measured in units of e 2 /h and L the system 
size, in the presence of Coulomb interactions for the case 
when a xx S> 1. The one-loop renormalization leads to 
the following equation 

f3{a xx ) = ^-l + {N 2 -l)T, (42) 

where N represents the degeneracy, which is 4 in 
graphene and 1 in TI. The term N/2 represents weak 
antilocalization, while the second term arises from the 
Coulomb interaction in the singlet channel and evidently 
suppresses the conductivity. It is independent of N since 
all the flavors are involved in screening. The last term 
is due to the Coulomb interaction in the multiplct chan- 
nel, with T an interaction parameter. In TI the last term 
vanishes and the scaling function turns out to be negative 
in the presence of electron-electron interactions, imply- 
ing a tendency towards localization. Since at small a xx 
the scaling function is positive due to topological protec- 
tion, it must cross zero at some critical point. The au- 
thors thus identify an interaction- induced critical point, 
concluding that interactions can fully localize topological 
surface states in 3D TI. 



VI. TRANSPORT EXPERIMENTS 

Experimental work on transport in topological insula- 
tors has begun to advance at a brisk pace, and is without 
doubt entering its heyday, in the way ARPES and STM 
work did two years ago (for a recent experimental review 
see Ref. 11070 It will help to summarize in the following 
what one expects to see in a sample in which bulk carriers 
have been eliminated and transport is due to the chiral 
surface states. When bulk carriers are still present one 
expects variable range hopping characterized by an ac- 
tivation energy. The semi-logarithmic plot of resistivity 
versus temperature in this case displays a clear plateau. 
Observation of a graphene-like minimum conductivity re- 
quires tuning through the Dirac point. Measuring the 
spin polarization of a current is not an easy experiment. 
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In contrast, weak localization and antilocalization both 
disappear in a weak magnetic field. A magnetic field also 
allows observations of Shubnikov-de Haas oscillations, 
which may be grouped into oscillations due to bulk states 
and oscillations due to surface states. The surface SdH 
oscillations are expected to vary as A(l/B) — e/(hn), 
where we recall that n is the 2D carrier density, and to 
show angular dependence: for surface transport one an- 
ticipates that SdH oscillations will only be present for 
B perpendicular to the surface. On the other hand, for 
bulk transport SdH oscillations are expected for all ori- 
entations of the magnetic field. The cyclotron mass, dis- 
cussed above under magnetotransport theory, is typically 
extracted from cyclotron resonance experiments. It is 
given by Ep/v F , where vp is the (constant) Fermi veloc- 
ity given by A/h. These latter considerations show that 
certain signatures of chiral surface states may be easier 
to observe using magnetotransport. Consequently, thus 
far significantly more experimental work has focused on 
magnetotransport than on non-magnetic transport. 

In the initial experimental efforts, it appeared impos- 
sible to identify any signatures whatsoever of the elusive 
surface states. One early avenue towards the removal of 
bulk carriers inv olve d compensation by doping with Ca. 
Checkelsky et a/.H2SI performed transport experiments on 
Ca-doped samples noting an increase in the resistivity, 
but concluded that surface state conduction alone could 
not be responsible for the smallness of the resistivities ob- 
served at low temperature. A subsequent paper by the 
same group^Sl claims to have observed surface conduc- 
tion, inferred from the sign change in the Hall coefficient. 
Such a sign change would indicates a change in the car- 
rier type from electrons to holes, and therefore the fact 
that the Dirac point has been crossed. The results are in- 
terpreted based on a model containing parallel bulk and 
surface conduction channels. The chemical potential is 
assumed to lie in the bulk gap, implying that bulk carriers 
are present due to thermal excitation across the gap. Cu- 
riously, the bulk gap is observed to be narrowed, greatly 
increasing the number of thermally excited bulk carriers. 
At the same time, the resistivity displays a power law 
variation as a function of temperature, whereas the acti- 
vated temperature dependence typical of variable range 
hopping is not observed. 

Analytis et alW^ used SdH oscillations and ARPES to 
probe the Fermi surface of single crystals of Bi2Sc3. The 
Bohr radius and Thomas-Fermi screening length place 
this material near the boundary of the metal-insulator 
transition. The SdH and ARPES probes agreed quanti- 
tatively on measurements of the effective mass and bulk 
band dispersion, and, in samples with high electron den- 
sity (10 19 cm -3 ), they agreed on the position of Ef- Yet 
discrepancies emerge regarding the location of £p in sam- 
ples with reduced carrier densities. For samples with 
carrier densities down to 10 17 cm -3 ARPES identifies the 
Dirac dispersion of the surface states, whereas SdH oscil- 
lations reveal a 3D Fermi surface, pinning the Fermi level 
in the bulk conduction band. None of the data displayed 



any thickness dependence. 

Eto et a/P321 reported angular-dependent Shubnikov- 
de Haas (SdH) oscillations due only to the residual 
3D Fermi surface in Bi2Se3 with an electron density of 
5 x 10 18 cm~ 3 . A significant finding by these authors is 
the fact that the bulk SdH oscillations themselves are 
anisotropic. A general guideline emerges that the obser- 
vation of anisotropic SdH oscillations by itself does not 
constitute a demonstration of surface transport. 

Steinberg et aZpH reported signatures of ambipolar 
transport analogous to graphene in Bi2Se3 samples with 
thicknesses of up to 100 nm, using a top-gate with a 
high-dielectric constant insulator. Ambipolar transport 
would be consistent with tuning between electron and 
hole carriers. Transport results were explained in terms 
of parallel surface and bulk contributions, and the dom- 
inant bulk contribution was subtracted so as to obtain 
the surface contribution. Surprisingly, the surface effect 
(i.e. an apparent minimum conductivity that would indi- 
cate ambipolar transport) was observed only for top-gate 
voltage sweeps and was absent during back-gate sweeps, 
a fact that to date rem ains unexplained. Further studies 
by the same group^^on weak antilocalization in Bi2Se3 
thin films revealed that the weak antilocalization contri- 
bution to the conductivity reflects the presence of various 
conduction channels - bulk as well as the top and bottom 
surfaces. 

Experimental studies of tetradymite materials have to 
date reported Landau levels but no quantum Hall effect, 
no conclusive measurement exists of the 7r-Berry phase 
in Shubnikov-de Haas oscillations. Cheng et alW* claim 
to have observed Landau levels analogous to graphene in 
100 nm-thick Bi2Se3 films, with Landau level energies oc 
V nB, where n is the Landau level index. The fitted data 
indicate an expected spin and valley degeneracy of one, 
yet puzzling features such as the absence of Landau levels 
below the Fermi energy and enhanced oscillations near 
the conduction band edge have not been resolved, and 
some arbitrariness is involved in extracting the Landau 
level filling factor v. Hanaguri et al¥^ report similar 
findings in the same material, with analogous features in 
the Landau level spectrum. 

Butch et oiP^I carried out SdH, ordinary Hall effect 
and reflectivity measurements on Bi2Se3 and concluded 
that the observed signals could be ascribed solely to bulk 
states, even at 50mK and low Landau level filling factor. 
Nominally undoped single crystals of Bi2Se3 with car- 
rier densities of ps 10 16 cm" 3 and mobilities in excess of 
2 m 2 /(Vs) were studied. Oscillations due to the surface 
states were expected at a very different frequency from 
those due to bulk states. The absence of a significant sur- 
face contribution to bulk conduction demonstrates that 
even in very clean samples, the surface mobility is lower 
than that of the bulk, despite its topological protect ion 
A similar conclusion was reached by Sushkov et aZP^ 
based on far-infrared cyclotron resonance and the Fara- 
day effect in Bi 2 Se 3 in magnetic fields up to 8 T, as well 
as by Jenkins et al. 116 based on the Kerr effect in the 
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same material. 

Taskin et aZP^3 studied bulk and surface transport in 
the material Bi1.5Sbo.5Te1.7Se1. 3, calculating that 70% of 
the total conductance arises from the surface. Data fit- 
ting yielded a surface (/-factor of approximately 20. Many 
intriguing unexplained findings were reported, including 
the observation that the surface states change as a func- 
tion of time. The scientists report quantum oscillations 
from Dirac holes and Dirac electrons in the same sample 
at different points in time over one month. This seems 
to be supported by the fact that the low-temperature 
Hall coefficient changes sign with time in thin samples, 
while the high temperature Hall coefficient is essentially 
unchanged as a function of time. Furthermore, the SdH 
maxima and minima were found to depend only on the 
perpendicular magnetic field. 

Qu et aZ.E^D studied quantum oscillations in Bi 2 Te3, 
claiming to resolve a surface electron current consistent 
with a Fermi velocity vf ~ 4 x 10 5 ms _1 , in agreement 
with ARPES measurements. The authors analyzed a va- 
riety of samples, of which some were metallic and some 
had a high bulk resistivity. Selective cleaving resulted in 
non-metallic crystals. Though in Bi2Te3 ep lies in the 
bulk valence band, the authors were able to manufacture 
both n- and p-type samples, which suggests a degree of 
competition and compensation between the bulk valence 
and conduction bands. Data from the weak-magnetic 
field Hall effect yield a surface mobility of 9000-10,000 
and substantially higher than in the bulk. SdH oscilla- 
tion data is shown in Fig. [8} In non-metallic samples, the 
SdH oscillation period depends only on the perpendicu- 
lar magnetic field, with oscillations not resolved for mag- 
netic field tilt angles q in the range 65 deg < q < 90 deg, 
whereas oscillations in the metallic sample survive up to 
90 deg. The Hall signal changes sign abruptly as Ef ap- 
proaches the Dirac point. However, the SdH oscillations 
appear to be weak, and were found by subtracting two 
large numbers. Even though the surface mobility is cal- 
culated to be very large and to exceed the bulk, even in 
the most resistive sample, the bulk conductance (i.e. to- 
tal signal) exceeds the surface conductance by a factor of 
« 300. The authors explain that the bulk hole density 
and mobility are not easily determined in the nonmetal- 
lic crystals. They also point out certain unconventional 
findings, for example, the transverse magnetoresistance 
of several non-metallic samples increases linearly with the 
magnetic field, a feature not seen in metallic samples. 

Coexistence of bulk a nd 2 D Fermi surfaces was also 
observed by Taskin et a/JliSI in their study of angular- 
dependent magnetoresistance and SdH oscillations in 
Bi .9iSbo.o9- SdH oscillations fall into two categories. 
Oscillations observed at low magnetic fields are believed 
to originate from the surface states. A different pattern of 
SdH oscillations becomes prominent at higher magnetic 
fields, which at present remains unexplained. 

In a beautiful experiment, Analytis et o/P^ recently 
detected the topological surface states of Bi 2 Se3 in which 
Sb was partially substituted for Bi to reduce the bulk 
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FIG. 8: Shubnikov-de Haas oscillations seen by Qu et al\ — ^ 
in two samples (labeled Ql and Nl in the article.) Figures (A) 
and (B) show the derivative dp xx /dH of the diagonal resistiv- 
ity with respect to the magnetic field H versus the perpendic- 
ular component H± of the magnetic field at selected tilt angles 
q. (A) Sample Ql. (B) Sample Nl. In (A), the minima lie on 
the vertical dashed lines consistent with a 2D Fermi surface, 
whereas in (B), the minima shift systematically with q. (C) 
The field position of the n = 3 Landau level peak for sample 
Ql (red circles) varies with q as 1/cosq (black curve), con- 
sistent with a 2D Fermi surface. (D) The period for sample 
Nl (blue circles) increases by 35 deg as q — > 90 deg, deviat- 
ing strongly from 1 / cos q (black curve) , consistent with a 3D 
Fermi surface. This figure has been adapted from Ref. 11181 



carrier density to 10 16 cm~ 3 . Results from this paper are 
shown in Fig. [9] The bulk carrier density was systemat- 
ically reduced by Sb doping until the quantum limit was 
reached, that is, the point where a magnetic field can 
collapse the bulk carriers to their lowest Landau level. 
Beyond this field surface state could be seen clearly. The 
carrier density determined by the low field Hall effect 
matches the Fermi surface parameters measured via SdH 
oscillations. When the field is rotated away from the 
crystal trigonal axis, the SdH frequency changes accord- 
ing to the morphology of the Fermi surface. SdH oscil- 
lations depend only on the perpendicular magnetic field 
and the oscillatory features grow with increasing field. 
Dips that are periodic in the inverse of the magnetic 
field in both the Hall derivative and the magnetoresis- 
tance align at all angles, leading the authors to conclude 
that the plateau-like features in the Hall and minima in 
the SdH oscillatins originate from the topological surface 
states. On the other hand, curiously, the quantization of 
the conductance is not evident. The 2D Fermi surface in 
the samples considered is small enough that the 2D quan- 
tum limit could be reached at pulsed fields of up to 60T. 
In this limit an altered phase of the oscillations was ob- 
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served, as well as quantum oscillations corresponding to 
fractions of the Landau integers, suggesting that correla- 
tion effects play an important role in the physics. These 
latter observations have not been fully explained. 

Xiu et re port experimental evidence for the mod- 
ulation of topological surface states by using a gate volt- 
age to control quantum oscillations in back-gate Bi 2 Tc 3 
nanoribbons. The bulk resistance increases below 50K 
as bulk carriers are frozen out, but little change in the 
resistance is measured as the temperature is reduced be- 
low 20 K, and a clear plateau of resistance versus tem- 
perature was not seen. Surface conduction was inferred 
from a small variation in the resistance below 4 K. Un- 
der zero gate bias at 1.4 K SdH oscillations are not seen. 
The magnetoresistance is linear with superimposed uni- 
versal conductance fluctuations. Although Bp moves, vp 
extracted from SdH oscillations remains constant, con- 
sistent with 2D transport. Up to half of the conductance 
is attributed to the surface states. The authors find that 
a gate voltage yields a significant enhancement of sur- 
face conduction, and the surface mobility may reach 5800 
cm 2 /(Vs), much higher than previously recorded. At the 
same time Aharonov-Bohm oscillations with a period of 
h/2e are reported. 

Chen et aZP^I studied the low-magnetic field weak an- 
tilocalization magnetoconductivity of Bi 2 Se 3 thin films, 
in which carrier densities can be tuned over a wide range 
with a back gate. Weak localization magnetoconductiv- 
ity was fitted to the well-known Hikami-Larkin-Nagaoka 
formula. For the case in which the carriers can be di- 
vided into an electron layer at the top surface and a 
hole layer at the bottom, the magnetoconductivity devi- 
ates strongly from the single-component Hikami-Larkin- 
Nagaoka equation. The authors found that, as the elec- 
tron density is lowered, the magnetotransport behav- 
ior deviates from the single-component description, and 
strong evidence exists for independent conducting chan- 
nels associated with the top and bottom surfaces. This 
description is found to be valid for surface electron den- 
sities in the range 0.8 — 8.6 x 10 13 cm -2 in samples in 
which the Fermi energy lies in the conduction band. 

A mor e recent experimental breakthrough by Sacepe 
et a/jA22l investigated surface transport in thin films of 
Bi2Se3 of thickness ps lOnm. Results from this work are 
shown in Figs. [10] and [IT] The authors fabricated super- 
conducting junctions on thin Bi2Se3 single crystals. A 
gate electrode allowed them to vary the electronic den- 
sity. In the presence of a sufficiently large perpendicular 
magnetic field, sweeping the gate voltage enabled obser- 
vation of the filling of Landau levels of the Dirac fermions. 
The Landau levels evolve continuously from electron-like 
to hole-like. At zero magnetic field and with the elec- 
trodes are in the superconducting state, a supercurrent 
appears, whose magnitude can be gate tuned, and is min- 
imum in correspondence of the charge neutrality point 
determined from the Landau level filling. While these 
results offer grounds for hope, they have not been repro- 
duced by other experimental groups to date. 



Kim et aZf^H investigated transport and weak an- 
tilocalization in Bi2Se3, finding a strong thickness- 
dependences over the thickness range 3 nm - 170 fim. 
The mobility was found to increase linearly with thick- 
ness in the thin film regime and saturated in the thick 
limit, indicating dominant bulk transport. The weak 
localization measurements were consistent with surface 
transport over two decades of thickness, yet the phase 
coherence length increased monotonically with thick- 
ness. This latter observation hints at strong bulk-surface 
hybridization. The same group later reported giant, 
thickness-independent surface transport in Bi2Se3 thin 
films with a surface-to-bulk conductance ratio believe 
to be ps 470% and surface-to-bulk conductivity ratio of 
11,000%. The films ranged in thickness from 2 to 2750 
quintuple layers. The Hall coefficient was measured to 
be nonlinear in the magnetic field, in agreement with a 
two carrier model with different mobilities. Over the en- 
tire thickness range the inferred surface density was unal- 
tered. Interestingly, this surface density persists down to 
two quintuple layers, which, as the scientists themselves 
remark, is the lowest thickness at which top and bottom 
surfaces can still be defined. 

Quantum oscillations were also recently seen in 
Bi2Te2SeP^3 In this material a high resistivity was 
recorded, in excess of 1 ffcm as well as variable-range 
hopping features, accompanied by SdH oscillations at- 
tributed to the topological surface states. The authors 
inferred a peculiar modified Berry phase from the SdH 
data, which is not fully understood. The surface accounts 
for 6 % of the total conductivity, implying a significant 
number of bulk carriers, and the calculated carrier den- 
sity from SdH oscillations does not match the low field 
Hall carrier density. These findings are ascribed to an im- 
purity band, yet the presence of an impurity band is puz- 
zling in well-ordered chalcogenide layers. For deep lev- 
els the carrier density needed to form an impurity band 
greatly exceeds the calculated bulk den sity. A different 
group working on the same material^-' reported SdH os- 
cillations in a sample with a bulk resistivity of 6 i7cm, 
and a surface mobility of 2800 cm 2 /(Vs), two orders of 
magnitude larger than the bulk mobility. 

Cho et aZf^H grew ultrathin (approximately three 
quintuple layer) n-type field-effect transistors of Bi2Se3. 
The devices have a clear OFF state and display acti- 
vated temperature-dependent conductance and energy 
barriers up to 250 meV, identifying them as conven- 
tional insu lator s. In a more recent breakthrough, the 
same group^H showed that Bi 2 Se3 surfaces in samples 
with thicknesses of < lOnm, with carrier densities of 
< 10 17 cm~ 3 are strongly electrostatically coupled. A 
gate electrode entirely removes bulk carriers and takes 
both surfaces through the Dirac point simultaneously. 
Ambipolar transport was observed with with well-defined 
p and n regions, together with a minimum conductiv- 
ity of 3 — 5e 2 //i per surface, reflecting the presence of 
electron and hole puddles. Their results are shown in 



Fig. 12 This experiment also measured the weak lo- 
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FIG. 9: Shubnikov-de Haas oscillations seen by Analytis et a/P^ Longitudinal resistance i?^ (a) and Hall resistance R xy (b) 
data traces as a function of magnetic field for indicated angles for a single crystal with carrier density n w 4 x 10 16 cm -3 (labeled 
sample El in the paper.) The i?^ traces are offset for clarity. The deviation from linearity of the low-field R xy indicates the 
3D quantum limit (B sa AT), which does not vary with angle. Beyond this limit, additional features in both R xy and R xx 
move smoothly up in field as the tilt angle 8 is increased. 6 = is defined to be when the field is perpendicular to the surface 
(parallel to the trigonal c-axis). (c) A smooth second-third order polynomial can be fitted as background and subtracted from 
the raw R xx , revealing oscillations that grow with field, (d) -d 2 R xx /dB 2 and (e) dR xy /dB (e) as a function of 1/B±, where 
B± = BcosO aligns all features associated with the 2D surface state. The vertical lines in c-e indicate the first three Landau 
levels n = 1, 2, 3 of the 2D state. This figure has been adapted from Ref. 120 
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FIG. 10: Maximum resistivity, a signature of ambipolar 
transport. The figure shows the dependence of the normal 
state resistance on the gate voltage V g . The maximum at 



-10V, provides a first indication of ambipolar transport. 



The V a dependence originates from the modulation of the con- 
ductivity of the surface close to the gate, S<j(V g ). The inset 
shows the calculated 8a (V g ). This figure has been adapted 
from Ref. fT2~3l 



calization magnctoresistance consistent with topological 
surface transport. 

Two recent papers from the Wiirzburg group have re- 
ported exciting developments in HgTe, a semimetal. The 



first was a study of zero-gap HgTe quantum wells 
which realizes a two-dimensional electron system with 
a single spin-degenerate Dirac valley. The quantum well 
was grown at a critical thickness where the bandgap van- 
ishes. Application of a perpendicular magnetic field re- 
sulted in the observation of the quantum Hall effect in 
this novel material. The quantum Hall plateaux dis- 
played the anomalous sequence characteristic of a (single- 
valley) Dirac systems. The second paper concerned the 
observation of the quantum Hall effect due to the topo- 
logic al s urface states of a strained 70 nm-thick HgTe 
layer This marks the first time that a well-defined 
Hall conductance from the surface states is observed ex- 
perimentally. The HgTe layer was strained by epitaxial 
growth on a CdTe substrate, which opens a gap in the 
bulk material. Strained bulk HgTe thus becomes a 3D 
TI. Contributions from residual bulk carriers to the trans- 
port properties of the gapped HgTe layer are negligible 
at mK temperatures. As a result, the sample exhibits a 
quantized Hall effect that results from the 2D single cone 



Dirac-like topological surface states, shown in Fig. 13 



In a related paper from the same group, Tkachov et 
aZP^l presented a theoretical and experimental study of 
the mobility of 5-12 nm-wide HgTe quantum wells with a 
finite gap. Whereas semiconductor heterostructures ex- 
hibit an increase in mobility with carrier density, a mo- 
bility maximum is seen in HgTe quantum wells. This is 
attributed to fluctuations in the well width, which lead 
to fluctuations in the mass of the carriers, a phenomenon 
that allows backscattering between states of opposite mo- 
menta. 

Chen et aZP^I investigated weak antilocalization in 
Bi2Se3 thin films epitaxially grown on SrTi03 substrates. 
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FIG. 11: 2D Landau levels for massless Rashba-Dirac 
fermions. The main panel shows the Landau level index N 
versus 1/B plotted for selected values of gate voltage, cor- 
responding to different concentration of holes (red symbols) 
and electrons (blue symbols). According to the relation char- 
acteristic for massless Dirac fermions, iV = 2g — | , where n 
is the density of charge carriers (electrons or holes), whereas 
for normal fermions N = jk. Examining the way N ex- 
trapolates for 1/B — > allows one to discriminate between 
massless Rashab-Dirac and normal fermions. The continuous 
lines are linear least-square fits to the data (N — a/B + b). 
They show that as 1/B — > 0, N extrapolates to close to —1/2 
as expected for massless Rashba-Dirac fermions (and not for 
normal fermions) . The inset shows the extrapolation value (b) 
obtained by least square fitting for approximately 20 differ- 
ent values of gate voltage (in the region with V g between —40 
and +15V, the oscillations are not sufficiently pronounced to 
perform the analysis.) The data show that b fluctuates close 
to —1/2 and not around 0, and, consistently, the average of 
b over gate voltage is (6) = —0.4 ±0.1, compatible only with 
massless Rashba-Dirac fermions rather than normal fermions. 
This figure has been adapted from Ref. 11231 



The dielectric constant of this substrate is extremely 
large (10000), enabling significant tunability of the car- 
rier density with a back gate. With the exception of 
very low electron densities, the low field weak antilocal- 
ization magnetoconductivity, shown in Fig. |14[ shows a 
very weak dependence on the gate voltage. These results 
indicate a heavily suppressed bulk conductivity at large 
negative gate voltages and a possible role of surface states 
in weak antilocalization. Even though a feature similar 
to the expected minimum conductivity was observed, a 
residual bulk conductivity cannot be ruled out based on 
present data. 

Weak antilocalization was also studied in thin films of 
Bi 2 Te3p22 The observed signal is believed to arise from 
the topological surface states due to the anisotropy of 
its dependence on the magnetic field. The weak antilo- 
calization signal is robust against deposition of nonmag- 
netic Au impurities on the surface of the thin films, but 
it is quenched by the deposition of magnetic Fe impuri- 
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FIG. 12: Minimum conductivity of a Bi2Se3 thin film. The 
figure shows the conductivity a xx and Hall carrier density 
measured in an AI2O3 capped device at T = 4.2K as a func- 
tion of back gate voltage. This figure has been adapted from 

Ref. mo 
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FIG. 13: 2D Landau levels in strained HgTe. The figure 
shows transport data on a strain-gapped 70-nm-thick HgTe 
sample. The Hall conductivity measured at 50 mK shows 
plateaus at the quantized values. The inset shows the Hall 
resistance R xy , together with the longitudinal resistance R xx . 
This figure has been adapted from Ref. 11311 



ties which destroy the it Berry phase of the topological 
surface states. Their results are shown in Fig. 
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Liu et aiP^ studied ultrathin Bi2Se3 films with thick- 
ness from one quintuple layer to six quintuple layers 
grown on sapphire. At low temperatures, the film re- 
sistance increases logarithmically with decreasing tem- 
perature, which is believed to represent an insulating 
ground state. The resistance increase as a function of 
temperature becomes more pronounced in thinner films. 
The sharp increase of resistance with magnetic field is in- 
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FIG. 14: Weak antilocalization correction to the conductiv- 
ity. The figure shows the magnetoconductivity Sa™^ 9 of Eq. 
( 36 1 as a function of magnetic field measured in various sam- 
ples. Figure (a) is for the sample labeled A in the article, and 
(b) and (c) are for the sample labeled B at different sets of 
gate voltages. Figures (d)-(f) are fitted curves of the magne- 
toconductivity of sample B at different gate voltages using the 
Hikami-Larkin-Nagaosa formula, given in the same equation. 
This figure has been adapted from Ref. 11331 



terpreted as weak antilocalization due to the topological 
surface states. The authors claim the insulating ground 
states is induced by electron-electron interactions, yet it 
is puzzling that electron-electron interactions can be ob- 
served when the dielectric permittivity of Bi 2 Se3 is ap- 
proximately 1 00. 

Wang et a/P^ studied thin films of Bi2Se3 with and 
without Pb doping. ARPES clearly identifies the topo- 
logical surface states in the Pb doped films, with ep in 
the bulk gap. The thickness of the films is smaller than 
the inelastic mean free path. Transport data show weak 
localization behavior, but the authors claim the tempera- 
ture and magnetic field dependences of resistance cannot 
be explained without appealing to electron-electron in- 
teractions. On the other hand, both the bulk and the 
surface participate in transport, and a clear separation is 
not possible in the samples investigated. 



VII. EXPERIMENTAL CHALLENGES 

The biggest obstacle to the observation of Z2 topolog- 
ical order in transport has been a factor that may be 
termed the bad metal syndrome. The existence of chiral 
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FIG. 15: Weak antilocalization correction to the conductiv- 
ity. The figure shows the normalized magnetoresistance of a 5 
nm Bi2Te3 film (a) and a 50 nm Bi2Te3 film (c) measured in 
tilted magnetic fields at T — 2K. The solid curve in (a) is a 
parabolic fit to the magnetoresistance data measured at zero 
tilt angle. The magnetoconductance as a function of the nor- 
mal magnetic field component with the magnetoconductance 
at zero tilt subtracted for the 5 nm sample (b) and for the 50 
nm sample (d). This figure has been adapted from Ref. 11341 



surface states is well established experimentally through 
angle-resolved photoemission (ARPES) and scanning- 
tunneling microscopy (STM) studies, yet in transport 
their contribution is overwhelmed by the bulk, which, due 
to unintentional defect-induced doping, is a low-density 
metal. To date no experimental group has produced 
a true topological insulator. Consequently, all existing 
topological insulator systems are in practice bulk-doped 
materials. 

One key aim of current TI efforts is the identification of 
a definitive signature of surface transport. Several recent 
experimental works point to a possible way forward, at 
least in the detection of surface transport, if not a perma- 
nent solution to the bul k dop ing problem. In recent gated 
thin film measurements f^H the two notable features are 
gate voltage tuning well past the Dirac point, and the 
resemblance of the conduction to that of dirty graphene 
dominated by charge puddles, signifying the removal of 
bulk electrons. Cho et al. also r ecen tly reported com- 
pletely insulating Bi2Se3 devices!^!! These researchers 
used mechanical exfoliations to fabricate ultrathin field- 
effect transistors of approximately three quintuple layers 
of Bi2Se3 on a Si/Si02 substrate. The devices are n-type 
and have a clear off state at negative gate voltage. The 
conductance shows an activated temperature dependence 
with a large activation energy up to 250 meV. 

Kong et a/P^ reported the synthesis and character- 
izations of ultrathin Bi2Te3 and Bi2Se3 few-layer tun- 
able nanoplates with thickness down to 3 quintuple lay- 
ers, using a catalyst-free vapor-solid (VS) growth mech- 
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anism. Optical images reveal thickness-dependent color 
and contrast for nanoplates grown on Si/Si02- Ultrathin 
TI nanoplates have an extremely large surface-to- volume 
ratio and can be electrically gated more effectively than 
the bulk form, potentially enhancing surface state effects 
in transport measurements. Low-temperature transport 
measurements of a single nanoplate device, with a high-k 
dielectric top gate, show decrease in carrier concentration 
by several times and large tuning of chemical potential. 

Teweldcbrhan et aZp^Sl obtained atomically thin crys- 
talline films of Bi 2 Te3 on Si/Si02 by exfoliation, which 
were characterized by atomic force microscopy and 
micro-Raman spectroscopy. The crystals have low ther- 
mal conductivity and high electrical conductivity. Re- 
markably, the films are non-stoichiomctric: altering the 
thickness and sequence of atomic planes enables one to 
change the composition and doping of the films. 

The above presentation demonstrates that, thus far, 
experimentalists have only succeeded in eliminating bulk 
carriers temporarily using a gate. The undesirable resid- 
ual carrier density can be reduced by doping, since the 
presence of dopants drains the excess carriers and brings 
the chemical potential down into the bulk gap. One po- 
tential way forward was opened by Hor et al.p^ who 
demonstrated that Ca doping (sa 1% substituting for Bi) 
brings the Fermi energy into the valence band, followed 
by the demonstration by Hsieh et alW^ that Ca doping 
can be used to tune the Fermi energy so that it lies in 
the gap. However, the Ca-doped material displays irreg- 
ular magnetotransport features that have not been ex- 
plained yet. Furthermore, a larger dopant concentration 
translates into a larger impurity concentration, and in 
the current TI materials this may cause the formation 
of impurity bands. A large impurity concentration may 
also lead to Anderson localization in the bulk, which is 
dangerous for TI. The bulk conductivity vanishes at zero 
temperature for an Anderson insulator, but bulk hop- 
ping contributes to transport at finite temperatures and 
may still dominate over the contribution from the sur- 
face states. Mechanically exfoliated thin TI films on a 
substrate have turned out to have large carrier densi- 
ties, hinting that further doping is introduced during the 
cleaving process or via interaction with the substrate. 

An additional experimental problem is the observa- 
tion th at t he atmosphere leads to rt-type doping of the 
surfaceP^Sllt is known that Bi2Sc3 acquires additional -re- 
type bulk doping after exposure to atmosphere, thereby 
reducing the relative contribution of surface states to the 
tot al o bserved conductivity. Kong et al. demonstrated 
thiJS^ via transport measurements on Bi2Se3 nanorib- 
bons. Systematic surface composition analyses by X- 
ray photoelectron spectroscopy reveal fast formation and 
continuous growth of native oxide on Bi 2 Se3 under ambi- 
ent conditions. In addition to n-type doping at the sur- 
face, such surface oxidation is likely the material origin of 
the degradation of topological surface states. Appropri- 
ate surface passivation or encapsulation may be required 
to probe topological surface states of Bi2Se3 by transport 



measurements. In this context, I will note that ARPES 
results obtained by different groups on the same mate- 
rial tend to show differences in the measured dispersion 
relations, with slightly different values of vp, and a sig- 
nificant time-dependence. 

Considerable work has focused on the detection of sur- 
face states using quantum oscillations. One recurring 
problem in current experimental reports of coexisting 
bulk and surface magnetotransport it the fact that in- 
terpretations typically rely on the angle dependence of 
the quantum oscillations, which in many cases is in fact 
very weak, and, until convincing proof is devised, may 
still reflect the anisotropy of 3D bulk transport. More- 
over, substantial disagreements exist between the Hall 
coefficient measured in different ways, and in many ex- 
perimental studies, although signatures of 2D electrons 
were reported, they were frequently accompanied by a 
host of unconventional extraneous findings, such as time- 
dependent Hall coefficients. Despite the reported obser- 
vation of SdH oscillations by many experimental groups, 
one salient issue, pointed out above in the discussion of 
the Hall effect is that the TI surfaces are all connected. 
Therefore in principle one does not expect the period of 
SdH oscillations to diverge. A full understanding of the 
surface SdH oscillations without doubt requires further 
investigation. 

These obstacles are exacerbated by the fact that, 
whenever an indication of surface mobility can be ob- 
tained, it is characteristically very low, often below 1000 
cm 2 /Vs. Below some threshold, when exfoliating bulk 
crystals, both the effective 3D carrier density and scat- 
tering go up, which is the primary reason why a sub- 
stantial effort has been devoted to the gating process. 
A working model that has gained popularity nowadays is 
that the cleaving process mechanically introduces various 
kinds of defects into the samples, an observation borne 
out by the fact that MBE sample growers have experi- 
enced great difficulty in obtaining carrier densities lower 
than 10 18 cm -3 . The band structure may be distorted 
near the surface due to space-charge accumulation. For a 
host of practical reasons, the procedure of removing the 
bulk electrons by reducing the sample thickness appears 
to bring in unwanted problems of its own. On a more 
fundamental level, surface electrons experience the ran- 
dom potential due to impurities, inhomogeneities, and 
roughness found at the interfaces. Topological insulator 
surfaces remain poorly understood, and experimentally 
the next big challenge may well turn out to be the task 
of engineering cleaner and more tractable surfaces. 

From the point of view of data interpretation, bearing 
in mind that the bulk conductance tends to overwhelm 
the surface conductance, subtracting the bulk from the 
total in order to obtain the surface contribution is tricky 
and may be unreliable. A small number that is obtained 
by subtracting two large numbers can often lead to er- 
roneous conclusions. In addition, many interpretations 
have hinged on a transport model that considers two 
types of carriers, which are taken to be surface and bulk 
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carriers, although the model itself does not distinguish 
2D from 3D transport. 

One final factor to bear in mind, though not a problem 
in itself, is that the Rashba-Dirac cone can be buried in 
the bulk conduction or valence bands. Its location within 
the bulk gap is not topologically protected. 

VIII. CONCLUSIONS AND OUTLOOK 

A feature of this review that should be plainly obvious 
is that the work covered spans a period of only three years 
- in fact, the bulk of it was reported within the last two 
years alone. The fact that a transport review is neces- 
sary after such a short period of time illustrates both the 
overwhelming fundamental interest in this field and its 
phenomenal pace of growth. Three-dimensional topolog- 
ical insulators have grown from non-existence to a vastly 
developed mature field involving hundreds of researchers 
practically overnight. Within this time span, chiral sur- 
face states started out as a mere theoretical concept, were 
predicted to exist in several materials, were subsequently 
imaged and have finally been seen in transport. The rise 
of topological insulators is following a close parallel to 
the rise of graphene a short time ago. Credit for this 
unprecedented development must be attributed to the 
tireless efforts of many experimental groups and to the 
state-of-the-art technology available in the present day. 

As the previous section made abundantly clear, sev- 
eral stumbling blocks remain in transport, particularly 
in experiment. The growth of new materials is a nontriv- 
ial issue, and obtaining high-quality samples where only 
the surface electrons can be accessed in transport has 
proved to be a difficult task. It is especially important 
to bear in mind that future work may initially be ham- 
pered by factors such as roughness and dirt inherent in 
solid-state interfaces. The theoretical results described 
in this review will become useful in discerning transport 
by the topological surface states only after a real topo- 
logical insulator material, that is, a system which is a 
bulk insulator, becomes available. It remains true that 
all current TI materials are effectively bulk metals be- 
cause of their large unintentional doping - at present, 
bulk carriers are only removed temporarily by gating. 
Such systems are, by definition, unsuitable for studying 
surface transport properties since bulk transport over- 
whelms any signatures of surface transport. Discussing 
TI surface transport in such bulk-doped TI materials re- 
tains some ambiguity, since it necessarily involves com- 
plex data fitting and a series of assumptions required by 



the necessity of distinguishing bulk versus surface trans- 
port contributions. Although such analyses are being 
carried out in bulk doped TI materials, in order to un- 
derstand the transport behavior, real progress is expected 
when surface TI transport can be carried out unambigu- 
ously, without any complications arising from the (more 
dominant) bulk transport channel. 

Coexisting bulk and surface transport remains a rel- 
evant problem to TI experiments at present, and will 
be until scientists can produce a material that does not 
require gating in order to shift the chemical potential 
into the bulk gap. It could also be useful in the case of 
very small gap TIs that have large numbers of thermally 
excited bulk carriers. Analysis of current experimental 
data often uses a simple one carrier model, rather than 
relying on a systematic two-carrier analysis as a function 
of thickness. Also, it is not clear whether minimalistic 
heuristic models based on parallel surface and bulk con- 
duction are valid in systems in which these two types of 
carriers are present. 

Undoubtedly, the next major tasks facing experimen- 
talists are getting the chemical potential in the gap with- 
out the aid of a gate, further experimental studies con- 
firming the existence of a minimum conductivity, the 
measurement of a spin-polarized current, the search for 
the anomalous Hall effect (whether half-quantized or not 
quantized) and the associated giant Kerr effect, and po- 
tentially measurement of the giant spin-Hall effect. Theo- 
retically, areas that have received little attention include 
a comprehensive determination of the role of electron- 
electron interactions, the Kondo effect and the associ- 
ated resistance minimum, the possibility of large extrin- 
sic spin-orbit interactions, a clear resolution of the weak 
localization problem and its extension to cover universal 
conductance fluctuations. 

The observations above should not be considered a pes- 
simistic assessment, since topological insulators are truly 
new materials, and materials science must explore ways 
to handle novel problems, such as accessing the surface 
states. That is necessarily a challenging undertaking. Al- 
though fine control over carrier densities remains difficult 
with existing bulk-doped material, considering that the 
field has only been active for three years, the progress reg- 
istered is astounding, and tempering an occasional impa- 
tience for instant success can only be fruiftul in the long 
run. 

I am extremely indebted to N. P. Butch, S. Das Sarma, 
Yongqing Li, Changgan Zeng, Zhong Fang, Junren Shi, 
Di Xiao, and Shun-Qing Shcn. This work was in part 
supported by the Chinese Academy of Sciences. 
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